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ABSTRACT 


>This  report  describes  the  theory  and  operation  of  a  computer  program 
referred  to  as  SPADE  (Signal  Processing  and  Data  Enhancement).  SPADE  is  a 
computer  program  for  the  numerical  processing  of  transducer  signal  data  available 
in  the  form  of  digital  computer  files  from  air-blast  structural  response  trials. 
In  addition,  the  program  can  graphically  compare  the  processed  data  with 
predictions  made  using  the  finite  element  structural  analysis  computer  programs 
VAST  and  ADINA. 


SPADE  is  a  fully  interactive  program  that  can  be  run  from  a  computer 
terminal.  The  Tektronix  TCS  PLOT  10  and  AGII  packages  are  used  to  optionally 
give  an  interactive  plotting  capability.  The  program  has  a  reasonably 
comprehensive  set  of  options  which  makes  it  easy  and  convenient  to  process  the 
voluminous  amount  of  data  that  can  be  generated  during  a  shock  wave-structure 
interaction  experiment.  With  reasonable  ease  the  program  can  be  adapted  to  the 
processing  of  any  type  of  transducer  output.  However,  the  application  of  SPADE 
described  in  this  manual  is  to  the  processing  of  data  from  pressure  transducers, 
accelerometers,  and  strain  gauges  on  simple  plates  and  stiffened  panels  exposed 
to  air-blast  shock  wave  loading.^ 
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RALSTON,  ALBERTA 

SUFFIELD  MEMORANDUM  NO.  1123 

SPADE  -  An  Interactive  Computer  Program  for  Signal  processing  ^nd  Data 
Enhancement  of  Air-Blast  Structural  Response  Measurements  • 

by 

R.  Houlston 
and 

J.E.  Slater 

1.  INTRODUCTION 

The  Defence  Research  Establ ishment  Suf field  (DRES)  is  currently  conducting 
research  on  the  effects  of  air  blast  waves  on  warship  structures  such  as 
stiffened  panels.  The  work  includes  extensive  testing  of  both  model  and  full 
scale  structures.  Since  these  experiments  generate  voluminous  amounts  of  data, 
an  efficient  data  processing  procedure  is  required.  This  report  describes  a 
computer  program  SPADE  (pignal  processing  £nd  pata  Enhancement) ,  which  has  been 
developed  for  the  numerical  processing  of  transducer  signal  data  available  in  the 
form  of  digital  computer  files  from  air  blast  structural  response  trials.  The 
program  can  also  process  data  files  created  with  the  finite  element  computer 
programs  VAST  [1],  ADINA  [2],  or  DRESDYN  [3],  for  subsequent  plotting  with 
associated  experimental  results  for  data  correlation  studies.  SPADE  is  fully 
interactive  and  has  been  written  with  the  intent  of  making  it  as  easy  as  possible 
to  process  the  voluminous  amount  of  data  that  can  be  generated  during  a  shock 
wave-structure  interaction  experiment. 

SPADE-PLOT  and  SPADE-PLOTR  [4]  are  also  fully  interactive  computer  programs 
for  graph  plotting  of  the  finite  element  and/or  experimental  data  files  created 
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with  SPADE.  The  Cal  comp  plot  routines  are  used  for  high  quality  graphs.  SPADE 
also  has  the  option  of  using  the  Tektronix  PLOT  10  TCS  and  AGII  plotting  routines 
for  interactive  graph  plotting  during  data  processing  if  a  graphics  terminal  is 
being  used. 

Background  information  and  the  capabilities  of  SPADE  are  described  fully  in 
Sections  2  to  6.  The  use  of  the  program  is  then  discussed  in  Section  7. 

A  series  of  sample  runs  which  illustrate  how  to  run  the  programs  on  a 
computer  terminal  is  given  in  a  separate  document  [5]. 


2.  DATA  ACQUISITION 

Figure  1  shows  a  block  diagram  to  illustrate  the  data  acquisition  procedure 
from  an  experiment  and  the  processing  involved.  The  target  structure  is  equipped 
with  transducers  such  as  pressure  gauges,  strain  gauges,  and  accelerometers. 

Data  from  these  transducers  are  recorded  by  the  Electronic  Design  and 
Instrumentation  Group  (EDIG)  on  magnetic  tape  in  analog  form.  In  addition,  data 
for  calibration  pulses  (which  will  be  discussed  subsequently)  is  recorded.  An 
EDIG  Data  Sheet  (EDS)  is  also  compiled  for  recording  information  relevant  to  the 
data  acquisition  and  analog  signal  conditioning.  The  EDS  is  discussed  in 
Appendix  A. 

On  completion  of  the  experiment,  the  next  step  is  to  digitize  the  analog 
signals.  This  process  is  handled  by  the  Computer  Group  (CG)  with  information 
provided  by  EDIG  on  the  EOS.  As  a  result  of  the  digitizing  process,  a  new 
magnetic  tape  is  created  containi..g  the  same  data  as  recorded  by  EDIG  but  in 
digital  form.  Both  Calibration  Pulse  Data  Files  (CPDF)  and  Signal  Data  Files 
(SDF)  are  stored  on  this  tape  and  are  given  definite  file  names  (Appendix  B2). 
Before  processing  by  SPADE  can  begin,  the  files  required  must  be  extracted  from 
this  tape.  SPADE  also  produces  other  files  during  processing  referred  to  as  tlie 
Transducer  List  File  (TLF)  and  the  Transducer  Database  File  (TDF)  which  will  be 
discussed  in  Section  6.1  and  Appendix  B. 
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In  addition  to  creating  the  digitized  data  files,  the  Computer  Group  will 
also  produce  a  special  data  sheet  referred  to  as  a  Computer  Group  Data  Sheet 
(CDS)  to  be  discussed  in  Section  5. 


3.  PRELIMINARY  DEFINITIONS 


3.1  Introduction 


The  structure  referred  to  in  Figure  1  could  be  of  any  type  but  the  main 
mte'-est  at  the  time  of  writing  is  in  simple  plates  or  stiffened  panels  that  are 
Subjected  to  an  air  blast  shock  wave.  In  order  to  define  a  few  terms  that  will 
be  used  in  describing  the  use  of  SPADE,  a  schematic  illustration  of  a  simple 
pi jte  exposed  to  a  blast  wave  and  instrumented  with  transducers  is  shown  in 


"he  Charge  is  suspended  above  the  geometric  centre  of  the  plate  such  that 

tht'  line  joining  the  centre  of  the  plate  and  the  centre  of  the  charge  is 

perpendicular  to  the  plate  surface.  The  length  of  this  line  is  referred  to  as 
the  STANDOFF.  By  definition,  the  SURFACE  of  the  plate  that  is  exposed  to  the 

blast  wave  is  referred  to  as  surface  T  to  denote  the  top  surface.  The  bottom 

surface  is  denoted  as  surface  B. 

3.2  Transducer  Codes 


In  general,  several  different  types  of  transducer  could  be  used  to  measure 
the  blast  wave  characteristics  and  the  associated  structural  response.  For  the 
plate  in  Figure  2  only  a  pressure  gauge,  an  accelerometer  and  a  strain  gauge  are 
shown.  There  are,  in  general,  several  transducers  of  the  same  type.  For  the 
purpose  of  transducer  identification,  each  transducer  is  given  a  specific 
TRANSDUCER  CODE  of  the  form  XN  where  X  is  a  single  letter  to  denote  the 
transducer  type  and  N  is  an  integer  usually  of  one  or  two  digits.  For  a  pressure 
gauge  X  is  P,  for  an  accelerometer  A,  and  for  a  strain  gauge  S.  N  will  usually 
be  taken  as  1  for  the  first  transducer  of  a  particular  type  and  will  then  be 
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incremented  by  1  for  each  of  the  other  transducers  of  this  type.  However,  N  can 
be  assigned  arbitrarily  subject  to  being  unique  for  each  transducer  of  a 
particular  type.  In  Figure  2  the  transducer  codes  are  P3,  A2,  and  SIS. 

3.3  Transducer  Coordinates 


In  addition  to  identifying  a  particular  transducer  uniquely  with  a 
transducer  code,  it  is  necessary  to  define  the  position  of  that  transducer  on  the 
structure.  This  is  done  with  the  TRANSDUCER  COORDINATES.  In  general  a 
rectangular  Cartesian  coordinate  system  can  be  used  to  define  the  spatial 
position  of  a  point  uniquely.  In  the  case  shown  in  Figure  2,  a  two  dimensional 
rectangular  Cartesian  coordinate  system  is  used  with  the  origin  at  the 
geometrical  centre  of  the  plate  and  the  axes  passing  through  the  middle  of  the 
plate  sides.  Specification  of  the  surface  (T  or  B)  and  the  coordinates  (XC,  YC) 
is  sufficient  to  define  the  position  of  a  pressure  gauge  or  accelerometer 
uniquely.  In  the  case  of  a  strain  gauge  an  additional  coordinate  referred  to  as 
ANGLE  is  required.  This  coordinate  is  the  angle  that  the  line  along  the  length 
of  the  gauge  makes  with  the  XC  axis  as  shown  in  Figure  3.  ANGLE  is  measured 
positive  as  shown  (counter-clockwise  from  the  direction  of  the  positive  XC  axis). 
It  is  sufficient  to  require  0  <  ANGLE  <  180“  since  no  direction  needs  to  be 
associated  with  the  line  M.  This  means  that  a  rotation  of  the  strain  gauge  by 
180"  so  that  A  replaces  B  and  B  replaced  A  in  Figure  3  does  not  change  the  strain 
gauge  coordinates  (XC,  YC,  ANGLE). 

3.4  Transducer  Channels 

Each  of  the  transducers  shown  in  Figure  2  produces  an  electrical  signal  of 
magnitude  proportional  to  the  particular  stimulus  that  it  is  designed  to  measure. 
In  the  case  of  a  strain  gauge,  for  example,  the  change  of  length  as  it  is 
stretched  with  the  structure  changes  the  electrical  resistance.  Thus,  the 
electrical  current  variation  with  time  in  a  circuit  that  contains  the  strain 
gauge  is  a  measure  of  the  stretch  of  the  strain  gauge  or  more  precisely  of  the 
strain  on  the  surface  to  which  the  gauge  is  attached.  The  electrical  signal  from 
each  transducer  is  transmitted  to  the  electronic  data  acquisition  system  along 
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wires  attached  to  the  transducer.  Each  wire  is  referred  to  as  a  TRANSDUCER 
CHANNEL  or  simply  CHANNEL.  Usually  a  transducer  will  only  be  connected  to  one 
channel  but  in  general  several  channels  could  be  connected  to  a  single 
transducer.  This  requirement  would  hold,  for  example,  if  the  signal  from  a 
specific  transducer  was  required  to  be  recorded  with  different  signal 
conditioning  requirements.  Several  simultaneous  records  would  then  be  made  for 
the  same  transducer. 

Each  channel  is  given  a  number.  In  Figure  2  each  transducer  has  only  one 
channel.  A2  is  connected  to  channel  7,  P3  to  channel  3  and  S15  to  channel  107. 
These  channel  numbers  are  arbitrary  but  within  the  following  restriction.  In  the 
EDIG  instrumentation  trailer  there  are  two  tape  recorders  which  can  be  either 
Honeywell  (H-101)  and/or  Ampex  (A-PR2200)  as  illustrated  in  Figure  2.  Each 
channel  going  to  a  Honeywell  recorder  must  be  assigned  a  unique  number  between  1 
and  99.  On  the  other  hand,  if  an  Ampex  recorder  is  used,  each  channel  must  be 
assigned  a  unique  number  from  101  to  199.  This  assignment  of  channel  numbers  is 
related  to  the  file  naming  convention  for  Raw  Data  Files  (RDF)  discussed  in 
Appendix  B2. 

3.5  Series,  Shot  and  Panel  Numbers 


In  order  to  identify  a  particular  experiment,  a  Series  Number  (SERIES)  and 
Shot  Number  (SHOT)  are  introduced.  These  numbers  must  identify  uniquely  an 
experiment  but  can  otherwise  be  arbitrary.  The  Series  Number  will  usually  be 
kept  the  same  for  a  collection  of  experiments  which  only  differ,  for  example,  in 
the  standoff  and  charge  weight.  However,  the  Series  Number  would  probably  be 
altered  if  the  test  facility  or  type  of  structure  was  changed.  The 
classification  into  Series,  however,  is  essentially  done  in  any  way  that  is 
convenient.  Experiments  within  each  Series  are  then  given  Shot  Numbers,  usually 
starting  at  1  and  then  being  incremented  by  1  for  each  test  in  the  Series. 

To  account  for  the  possibility  of  more  than  one  panel  being  involved  in  a 
given  Series  and  Shot,  a  Panel  Number  (PANEL)  is  introduced.  This  number  is 
assigned  to  a  specific  panel  to  identify  it.  It  is  possible  that  several  shots 
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could  be  associated  with  the  same  Panel  in  one  or  more  Series,  but  the 
combination  (SERIES,  SHOT,  PANEL)  will  provide  unique  identification. 


4.  TRANSDUCER  ANALOG  SIGNALS 

4.1  Introduction 

As  mentioned  in  Section  2  and  illustrated  in  Figure  1,  the  transducer  analog 
signals  recorded  by  EDIG  must  be  converted  to  digital  form  before  numerical 
processing  can  begin.  The  recording  and  digitizing  processes  will  now  be  briefly 
summarized  in  order  to  define  certain  terms  that  occur  on  the  EDS  and  CDS. 

Figure  4  illustrates  in  graphical  form  the  digitizing  process. 

4.2  Recording  of  Analog  Signals 

Each  of  the  Ampex  or  Honeywell  tapes  has  14  tracks.  Twelve  of  these  tracks 
can  be  used  to  record  actual  transducer  signals  as  illustrated  in  Figure  4. 

Either  one  transducer  signal  is  recorded  on  each  track  or  multiplexing  is  used 
whereby  up  to  five  transducer  signals  can  be  recorded  on  each  track.  The  other 
two  tracks  are  used  generally  to  store  timing  and  firing  reference  information. 

During  an  experiment  such  as  that  illustrated  in  Figure  2,  the  following 
sequence  of  events  occurs.  About  ten  seconds  prior  to  charge  detonation,  both 
tape  drives  are  started  for  recording  analog  calibration  pulse  and  transducer 
signals.  Approximately  3  seconds  (±100  ms)  prior  to  the  charge  detonation, 
rectangular  voltage  calibration  pulses  of  known  magnitudes  (lasting  about  1 
second)  are  placed  on  the  tracks  of  the  tape.  The  magnitudes  of  these  pulses 
will,  in  general,  be  different  for  each  transducer  and  are  preset  depending  upon 
the  transducer  characteristics.  Following  this  calibration  signal  recording,  the 
signals  from  the  various  transducers  are  then  recorded  on  the  same  tracks  of  the 
magnetic  tape  as  the  associated  calibration  pulses. 
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At  the  time  of  charge  detonation,  a  rectcngular  voltage  pulse  recorded  on 
track  13  is  started  and  lasts  for  approximately  1  second.  Due  to  the  time 
required  for  the  blast  wave  to  reach  the  transducers  from  the  charge,  there  will 
be  a  delay  of  a  few  milliseconds  from  the  start  of  the  Det  Zero  Pulse  on  track  13 
and  the  detection  of  a  signal  from  the  transducers. 

Track  14  on  the  analog  tape  is  reserved  for  the  Inter-Range  Instrumentation 
Group  (IRIG)  signal.  The  IRIG  signal  comes  from  the  GOES  satellite  and  allows  a 
real  time  base  to  be  established  for  all  transducer  signals.  Usually  the  B  type 
signal  is  used  which  allows  the  time  at  any  point  on  the  tape  to  be  determined  to 
within  1  ms. 

4.3  Digitizing  the  Analog  Signals 

The  digitizing  procedure  consists  of  scanning  the  analog  tapes  and 
converting  the  analog  signals  to  digital  form  (see  Figure  4).  As  the  tape  moves 
past  the  tape  head,  samples  from  a  selected  track  are  taken  at  a  known  rate. 

This  rate  can  be  pre-selected  and  is  usually  specified  as  the  number  of  points 
per  second.  For  purposes  of  reference  this  will  be  referred  to  as  DRATE. 

From  each  track  on  the  analog  tape,  two  files  will  usually  be  created,  one 
is  the  Calibration  Pulse  Data  File  (the  CPDF)  and  the  other  is  the  transducer 
Signal  Data  File  (the  SDF)  as  indicated  in  Figure  4.  The  names  given  to  these 
files  depends  upon  the  experiment  identification  numbers  (SERIES  and  SHOT),  the 
Channel  Number  (CHANNEL),  whether  it  is  a  CPDF  or  SDF,  and  the  type  of  tape  drive 
(Honeywell  or  Ampex)  used  to  record  the  analog  signals.  This  file  naming 
convention  for  raw  data  files  is  discussed  in  Appendix  B2. 

The  digitizing  process  can  be  suimnarized  in  the  following  series  of  steps: 

1.  Scan  the  analog  tape  for  the  EVENT  from  the  time  decoded  from  the  IRIG  B 

signal.  The  event  refers  to  the  known  time  of  charge  detonation. 
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2.  Back  up  the  tape  for  about  5  seconds  and  scan  for  the  start  of  the 
calibration  pulse  on  the  track  corresponding  to  the  transducer  of  interest. 

3.  Back  up  the  tape  again  to  allow  a  predetermined  amount  of  lead  time. 

4.  Start  digitizing  the  calibration  pulse  into  the  CPDF  at  approximately  1666 
points  per  second  (1.666  kHz)  and  record  the  real  time  at  the  start  of 
digitizing  as  determined  from  the  IRIG  B  signal  (time  tj  in  Figure  4). 

5.  Stop  digitizing  the  calibration  pulse  when  the  required  number  of  samples 

have  been  taken  (usually  5000).  Since  the  calibration  pulse  widths  are 

approximately  1  second  (Figure  4),  this  number  of  digitized  points  will  give 
approximately  1  second  of  baseline  in  front  of  and  behind  the  calibration 
pulse. 

6.  Scan  the  track  on  the  analog  tape  for  the  start  of  the  transducer  signal. 

This  usually  occurs  about  2  seconds  after  the  end  of  the  calibration  pulse. 

7.  Back  up  the  tape  to  allow  a  preset  amount  of  lead  time. 

8.  Start  digitizing  the  transducer  signal  at  DRATE  samples  per  second  into  the 
SDF  and  record  the  real  time  value  which  corresponds  to  the  point  where 
digitizing  is  started  (time  t^  in  Figure  4). 

9.  Stop  digitizing  the  transducer  signal  when  the  required  length  of  signal  and 
number  of  sample  points  have  been  digitized. 

During  the  digitizing  process,  the  voltage  amplitudes  taken  from  the  analog 
tape  are  converted  to  numbers  rounded  to  the  nearest  integer  value.  The  size  of 
number  that  a  one  volt  amplitude  from  the  analog  tape  will  give  depends  upon  the 
gain  setting  (GAIN).  If  GAIN  =  1,  then  1  Volt  will  produce  the  number  204.8 
(although  this  would  be  stored  as  205).  As  the  gain  is  increased,  however,  this 
number  will  increase  in  proportion.  Thus,  if  GAIN  =  2,  then  a  1  Volt  analog 
sample  would  produce  the  number  409.6  (stored  as  410). 
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5.  THE  COMPUTER  GROUP  DATA  SHEET  (CDS) 


The  CDS  contains  all  relevant  information  required  by  SPADE  plus  some 
additional  information  for  record  purposes.  Table  5.1  shows  a  sample  CDS  with 
information  that  is  used  for  the  sample  run  to  be  presented  in  Reference  5. 

There  is  one  CDS  for  each  experiment  (i.e..  Series,  Shot  and  Panel). 

The  sample  CDS  indicates  that  the  experiment  referred  to  is  SHOT  4  in 
SERIES  1  which  was  done  on  31  March  1983.  The  experiment  was  for  a  charge  of 
113  g  suspended  3  m  above  the  target  plate  (see  Figure  2)  which  is  identified  as 
PANEL  2.  The  EVENT  indicates  the  date  and  time  of  the  experiment  in  the  format 
D/HM;S.MS  where  D  is  an  integer  of  3  digits  indicating  the  day  of  the  year,  H  is 
an  integer  of  two  digits  giving  the  hour  of  the  day,  M  is  an  integer  of  two 
digits  giving  the  minutes  into  the  hour,  S  is  an  integer  of  two  digits  giving  the 
number  of  seconds  into  the  minute,  and  MS  is  an  integer  of  3  digits  giving  the 
number  of  milliseconds  into  the  second.  The  accuracy  of  the  event  is  given  to 
the  nearest  millisecond. 

The  various  columns  of  the  CDS  in  Table  5.1  provide  information  for  each 
transducer  used  in  the  experiment.  Columns  1  and  2  give  the  Channel  Numbers  and 
tranducer  Codes,  respectively.  Three  of  these  transducers  (P3,  A2,  and  S15)  are 
illustrated  in  Figure  2.  Columns  3,  4,  and  5  give  the  transducer  coordinates 
(XC,  YC,  ANGLE).  XC  and  YC  are  specified  in  millimetres  and  ANGLE  in  degrees. 
Column  6  gives  the  surface  (T  or  B).  Columns  7-12  give  the  Digitizing  Start 
Time  for  the  transducer  signals  (time  in  Figure  4).  For  PI,  for  example,  the 
start  time  of  digitizing  was:  Year  (Y)  =  1983,  Day  (D)  =  90th,  Hour  (H)  =  14, 
Minute  (M)  =  9,  Second  (S)  =  59,  Millisecond  (MS)  =  775.  Column  13  gives  the 
digitizing  rate  (DRATE)  for  the  SDF  and  is  required  by  SPADE.  Column  14  gives 
the  value  of  CALIBRATE.  This  is  the  height  in  real  units  of  the  Calibration 
Pulse  for  the  transducer  [usually  pounds  per  square  inch  (PSI)  for  pressure 
gauges,  microstrain  (^e)  for  strain  gauges,  and  gravities  (g)  for 
accelerometers  I . 
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Columns  15  and  16  in  Table  5.1  give  the  GAIN  and  the  Gauge  Calibration 
Constants,  respectively.  These  quantities  are  only  required  by  SPADE  in  the 
exceptional  case  that  a  CPDF  does  not  exist. 


6.  "SPADE"  PROCESSING 


6.0  Introduction 


In  Figure  1,  the  purpose  of  SPADE  in  the  processing  of  digitized  transducer 
signals  is  surntnarized.  Before  processing  can  begin,  the  required  CPDF  and  SOF 
must  be  available  on  computer  disc  files.  In  addition,  all  information  on  the 
COS  must  be  available. 

There  are  a  total  of  12  options  available  in  SPADE.  Within  restrictions 
that  will  be  pointed  out  in  subsequent  discussions  in  this  section,  all  options 
can  be  executed  independently.  Also,  the  execution  of  the  required  options  can 
be  done  in  a  single  run  of  SPADE  or  spread  over  several  runs.  The  program  is 
fully  interactive  and  has  been  designed  to  be  as  simple  and  convenient  to  use  as 
possible.  The  program  also  has  a  comprehensive  set  of  error  diagnostics. 

In  the  following  subsections,  an  overview  of  the  various  options  will  be 
given  together  with  any  required  theory.  The  actual  use  of  the  program  is 
discussed  more  fully  in  Section  7,  and  a  large  selection  of  sample  runs  is 
presented  in  Reference  5. 

SPADE  is  executed  on  the  Honeywell  DPS  8/70  by  giving  the  CP6  command 

! SPADE. A 


where  A  is  the  account  number  containing  the  program.  SPADE  will  then  request 
certain  information  including  the  SERIES,  SHOT,  and  PANEL  to  identify  the 
expe-'iment.  In  addition,  a  list  of  12  possible  options  will  appear.  These 
options  are: 
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Option  Description 

1  Set  up  the  TLF  and  the  TDF  for  a  new  SHOT. 

2  Determine  the  SIGNAL  SCALE  FACTOR  from  a  CPDF. 

3  Scale  the  transducer  SIGNAL  DATA  FILE. 

4  Define  t  =  0  from  pressure  signals. 

5  Determine  pressure  pulse  characteristics. 

6  Integrate  acceleration  to  obtain  relative  velocity. 

7  Integrate  velocity  to  obtain  relative  displacement. 

8  Obtain  acceleration,  velocity,  and  displacement  based  on  fit  to 
displacement. 

9  FFT  and/or  filtering  of  a  transducer  signal. 

10  Create  averaged  pressure  file  and  VAST  load  file. 

11  Extract  displacement,  velocity,  acceleration  from  VAST  file  for 
plotting. 

12  Horizontal  or  vertical  shift  with  optional  rescaling. 

Before  discussing  each  of  these  options  in  detail,  the  Transducer 
List  File  (TLF)  and  the  Transducer  Database  File  (TDF)  will  be  discussed  briefly. 
The  format  and  names  of  these  files  are  given  in  Appendix  B. 

The  TLF  is  used  for  listing  purposes  only  and  contains  the  Channel, 
Transducer  Code,  Transducer  Coordinates,  and  the  Surface  Identification  Code 
(SUR)  from  columns  1  to  6  of  the  CDS  as  defined  in  Sections  3  and  5.  The  TDF 
also  contains  these  quantities  as  well  as  the  Digitizing  Start  Date/Time,  DRATE 
and  CALIBRATE  from  columns  7  to  14  of  the  CDS  (see  Section  5).  Thus  all 
quantities  that  are  given  on  the  CDS  (except  the  GAIN  and  Calibration  Constant  in 
columns  15  and  16)  are  stored  on  the  TDF  during  execution  of  Option  1.  If  the 
CPDF  does  not  exist,  the  GAIN  and  Calibration  Constant  are  input  during  Option  1 
and  immediately  converted  to  the  scale  factor  (E)  to  be  used  to  scale  the  SDF  to 
engineering  units.  E  is  then  stored  on  the  TDF. 

In  addition  to  the  CDS  data  values  (which  are  known  before  the  start  of  any 
SPADE  processing)  there  are  certain  quantities  which  are  required  to  be  known  but 
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are  only  determined  after  certain  options  have  been  completed.  These  additional 
quantities  are  listed  in  Table  6.1  together  with  the  option  number  in  which  they 
are  determined.  They  will  be  described  in  the  following  sub-sections. 

Table  6.1 


Additional  Quantities  on  the  TUF  Determined  by  Various  Options. 


Quanti ty 

Option 

Scale  Factor  E. 

2  (If  CPDF  exists) 

1  (If  CPDF  does  not 

exi St) . 

Calibration  Pulse 

2 

Definition  Points  P^  to  P^. 

Signal  Definition 

3,  5,  8 

Points  Pj  to  P^. 

Point  Number  on  the  SDF 

4 

for  t  =  0  (Po). 

Decay  Constant  (b^)  and 

5 

maximum  pressure  (Q^)  for 

Pressure  Tranducer  Signals 

from  a  simple  exponential 

curve  fit. 

Maximum  Impulse  (IMAX)  for 

5 

a  Pressure  Transducer 

Signal  from  integration  of 

experimental  data. 
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6.1  Option  1:  Set  Up  Transducer  Database  Files  for  a  New  Shot 

This  option  must  be  executed  before  any  other  processing  by  SPADE  can 
commence.  The  main  purpose  is  to  create  the  TLF  and  TDF  files  for  the 
experiment  under  consideration.  Once  these  files  have  been  created,  the  CDS  will 
no  longer  be  required.  All  information  that  the  program  requires  for  further 
options  for  the  particular  experiment  is  stored  on  the  TDF.  As  other  options  are 
executed  this  file  will  be  updated.  Appendix  B,  Table  B4.1,  contains  a 
description  of  the  format  of  this  file.  Quantities  in  Table  6.1  that  are  updated 
by  other  options  are  written  as  zero  to  the  TDF  by  Option  1  (with  the  exception 
of  E  when  a  CPDF  does  not  exist  as  noted  in  Table  6.1).  Option  1  will  request 
the  following  data  for  each  transducer:  CHANNEL,  TCODE,  XC,  YC,  ANGLE  (if 
applicable),  SUR,  Y,  0,  H,  M,  S,  MS,  and  DRATE  defined  in  columns  1  to  13, 
respectively,  of  the  COS  (see  "Terminology"  at  the  start  of  this  manual  for  a 
review  of  the  meaning  of  these  terms).  The  remaining  quantities  that  Option  1 
will  request  are  given  in  columns  14  to  16.  There  are  two  possibilities  here 
depending  upon  whether  the  CPOF  exists.  Usually  it  will  exist  but  if  SPADE  is 
unable  to  locate  it  in  the  directory  specified,  the  user  will  be  asked  to  specify 
whether  it  exists  anywhere.  The  information  requested  by  SPADE  will  then  depend 
upon  whether  or  not  it  exists  and  can  be  accessed. 

If  the  CPDF  exists,  only  the  quantity  Calibrate  from  Column  14  will  be 
requested.  Nothing  from  Column  15  or  16  will  be  required.  The  Calibrate  value 
and  the  information  on  the  CPDF  will  then  allow  Option  2  (to  be  discussed  in  the 
next  subsection)  to  determine  the  Scale  Factor  (E)  required  to  convert  the  SUF  to 
real  units  (RU). 

If  the  CPDF  does  not  exist,  sufficient  information  must  be  given  to  the 
program  to  enable  E  to  be  determined  directly  during  execution  of  Option  1. 

This  is  done  by  supplying  the  GAIN  and  Calibration  Constant  (C)  from  columns  15 
and  16,  respectively.  How  E  is  obtained  from  these  values  will  now  be 
expl ained. 
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As  discussed  in  Section  4.3,  one  volt  on  the  analog  tape  will  produce 
204.8  X  GAIN  digitized  units  (DU)  on  the  CPDF  (see  Figure  4).  The  Calibration 
Constant  C,  however,  is  the  number  of  real  units  of  stimulus  that  must  be 
applied  to  a  gauge  to  give  one  volt  of  output.  For  a  strain  gauge,  C  could  be 
expressed  as  microstrain  per  volt  and  for  a  pressure  gauge  as  PSI  per  volt.  Thus 
for  1  volt,  C  real  units  (ftU)  will  give  204.8  x  GAIN  DU's.  The  Scale  Factor  E  is 
the  number  of  RU's  per  DU  and  hence 


204.8  X  GAIN 


(6.1.1) 


Specification  of  GAIN  and  C  from  columns  15  and  16,  respectively,  thus 
allows  E  to  be  determined  directly  by  SPADE  during  this  option.  E  is  then  stored 
on  the  TUF. 

The  final  point  that  is  worth  mentioning  with  Option  1  is  the  flexibility 
that  is  allowed  in  the  building  of  the  TLF  and  TDF.  In  practice  there  could  be 
several  hundred  transducers  each  of  which  will  require  data  to  be  entered  from 
the  COS.  For  this  reason,  special  features  have  been  built  into  the  program  for 
the  user's  convenience.  These  features  are: 

a)  Corresponding  identical  entries  need  not  be  repeated.  Only  the  terminal 
return  key  needs  to  be  depressed. 

b)  On  completion  of  entering  the  data  for  a  transducer,  the  program  can  be 
terminated.  When  the  user  decides  to  carry  on,  the  building  of  TLF  and  TDF 
can  continue  from  the  point  of  termination. 

c)  If  the  same  structure  and  transducers  are  to  be  used  in  several  experiments, 
the  TLF  and  TDF  do  not  have  to  be  created  in  full  for  each  experiment.  Once 
created  for  the  first  experiment,  new  TLF  and  TDF  can  be  created  for 
subsequent  experiments  using  information  in  the  previous  TLF  and  TDF. 

Entries  are  still  required  for  each  transducer,  but  quantities  such  as  the 
transducer  coordinates  are  not  required  to  be  re-specified.  The  restart 
option  described  in  b)  also  applies  in  this  case. 
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The  use  of  Option  1  is  fully  demonstrated  in  a  sample  run  in  Reference  5. 

6.2  Option  2:  Determine  the  Signal  Scale  Factor  from  a 
Transducer  Calibration  Pulse 

6.2.1  Introduction 


The  main  objective  of  Option  2  is  to  determine  the  Signal  Scale  Factor  (E) 
required  by  Option  3  to  convert  the  data  on  the  SDF  to  real  units.  As  mentioned 
in  Section  4.2  and  illustrated  in  Figure  4,  a  rectangular  voltage  pulse  of  known 
magnitude  is  placed  on  the  analog  tape  track  for  each  transducer.  There  are  a 
certain  number  of  real  units  of  actual  stimulus  that  the  gauge  must  experience  to 
generate  this  voltage.  This  number  of  real  units  is  derived  from  the 
manufacturer's  calibration  data  for  the  gauge  and  is  stored  in  column  14  of  the 
CDS.  This  quantity  will  be  referred  to  as  CAL. 

Once  the  calibration  pulse  has  been  digitized  as  discussed  in  Section  4.3, 
the  calibration  pulse  will  have  a  certain  height  in  digitized  units,  h(DU),  as 
shown  in  Figure  4.  h(DU)  is  equivalent  to  CAL  real  units.  Hence, 

CAL 

E  =  -  (6. 2. 1.1) 

h(DU) 

is  the  SDF  scale  factor  (the  number  of  real  units  RU  per  digitized  unit  DU). 

Since  the  same  recording  and  digitizing  procedure  is  used  for  the  transducer 
signal  on  the  same  track  as  the  calibration  pulse,  it  follows  that  E  can  be  used 
to  convert  the  digitized  transducer  signal  data  on  the  SDF  to  real  units.  The 
result  of  this  conversion  will  be  stored  on  a  file  referred  to  as  the  Scaled 
Signal  Data  File  (SSDF).  The  file  naming  convention  for  these  files  is  presented 
in  Appendix  B2. 

The  main  objective  of  Option  2  is  thus  to  calculate  h(DU)  from  the  CPDF  for 
a  transducer.  This  value  can  then  be  used  to  calculate  E  which  in  turn  is  used 
by  Option  3  to  convert  the  data  on  the  associated  SDF  to  real  units  with 
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background  correction.  The  results  of  processing  with  Option  3  are  then  stored 
on  the  SSDF. 

In  the  rare  case  in  which  the  CPDF  does  not  exist,  E  will  have  been 
determined  in  Option  1.  Executing  Option  2  in  this  case  will  result  in  an  echo 
of  the  value  of  E  stored  on  the  TDF.  The  sample  run  in  Section  2.2  of 
Reference  5  gives  a  full  illustration  of  the  use  of  Option  2.  Sections  6.2.2  to 
6.2.4  below  give  the  details  involved  in  the  calculation  of  h(DU). 

6.2.2  Definition  of  the  Calibration  Pulse  Definition  Points 

Figure  5  shows  a  schematic  illustration  of  a  digitized  calibration  pulse 
where  the  digitized  points  have  been  joined  by  straight  lines.  Typically  the 
total  number  of  digitized  points  on  the  CPDF  is  5000.  For  the  purposes  of 
reference,  SPADE  uses  a  set  of  six  Calibration  Pulse  Defintion  Points  (CPDP) 
denoted  as  P  ,  P  ,,  P^,  P^^,  P,^  and  P^^.  These  points  will  now  be  defined. 

Let  it  be  understood  that  by  Point  Number  (PN)  we  mean  the  actual  position 
number  of  a  point  on  a  digitized  file.  Thus  PN  1  (or  simply  point  1)  is  the 
first  point  on  the  file.  Similarly,  if  there  are  5000  points  on  file,  PN  5000 
will  be  the  last  point.  P^  to  P^  will  each  have  a  value  equal  to  a  particular 
Point  Number. 

Points  P;  to  P^,  inclusive,  will  be  defined  by  SPADE  as  a  set  of  points 
before  the  Calibration  Pulse  (CP)  that  are  to  be  used  for  background  correction. 

A  set  of  such  points  joined  by  straight  lines  is  illustrated  in  Figure  5. 

Usually  Pj  will  be  set  to  1  to  denote  the  first  point  on  the  file.  Point  must 
define  the  start  of  the  pulse  rise  as  closely  as  possible.  It  should  not 
identify  a  point  on  the  pulse  rise  such  as  X  in  Figure  5. 

Points  P^  to  P^,  inclusive,  will  be  defined  by  SPADE  as  the  set  of  points  on 
the  calibration  pulse  top.  P  should  identify  the  end  of  the  pulse  rise  and 
again  should  not  be  assigned  to  a  point  such  as  X  which  is  on  the  pulse  rise.  P^ 
will  usually  turn  out  to  be  just  1  greater  than  P^  since  the  rise  time  will 
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usually  be  less  than  the  time  difference  between  two  consecutive  points  given  by 
the  reciprocal  of  the  digitizing  rate.  Point  should  define  the  start  of  the 
pulse  drop. 

Points  Pg  to  P^  inclusive  will  be  defined  by  SPADE  as  a  set  of  points  after 
the  calibration  pulse  that  are  to  be  used  for  background  correction.  Usually  P^ 
will  be  set  to  the  last  point  on  the  file.  Point  P^  must  define  the  end  of  the 
pulse  drop  as  close  as  possible.  P^  will  also  usually  turn  out  to  be  1  greater 
than  P,  . 

6.2.3  Determination  of  the  Calibration  Pulse  Definition  Points 

SPADE  will  automatically  set  P^  to  the  first  point  on  file  and  P(^  to  the 
last  and  scan  the  CPDF  for  P,,  to  P^  .  The  method  used  to  determine  P.  to  P^  will 
now  be  described. 

Let  YMAX  be  the  maximum  value  recorded  on  the  CPDF  and  YMIN  the  minimum. 

Let  a  parameter  NOISE  be  defined  as  the  background  to  signal  ratio  in  percent. 
Then  the  noise  band  aB  is  defined  as 


aB  =  (NOISE/100)  X  (YMAX  -  YMINi 


(6. 2. 3.1) 


Let  Y.,  and  Y^  be  the  data  values  of  the  N'th  and  M'th  points  on  the  CPDF. 

N  n 

Then  Y^^  and  Yj^^  are  not  considered  to  be  significantly  different  if  the  following 
condition  holds. 


Yn  -  Y,  S  aB 


(6. 2.3. 2] 


Equation  (6. 2. 3. 2)  just  states  that  if  two  data  points  do  not  differ  in  value  by 
at  least  the  noise  band  then  the  difference  is  not  considered  to  be  significant. 

Point  is  determined  by  scanning  the  CPDF  from  point  1.  After  each  point 

“  N-1 

N  with  data  value  Y^  is  read,  it  is  compared  with  the  baseline  average  Y  j  of 
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~KJ_1 

the  Yj  values  for  points  I  =  1  to  N  -  1.  Since  Yj  will  be  approximately  the 
mean  line  through  the  noise  band,  if 

Y,,  -  Y  ^  (6. 2. 3. 3) 

NX  2 

then  point  N  is  considered  to  be  in  the  background  to  the  left  of  the  calibration 
pulse.  Otherwise  the  amplitude  of  point  N  is  considered  to  be  significantly 
~  N-1 

different  than  Y  ^  .  Suppose  Yj^  is  significantly  different.  In  this  case  is  it 

actually  on  the  pulse  rise  or  is  it  on  just  an  extraneous  spike  in  the  data? 

SPADE  checks  this  by  requiring  that  the  next  9  points  also  be  significantly  above 

the  average  Y  If  this  auxiliary  condition  is  not  satisfied  then  point  N  is 

not  considered  to  be  on  the  pulse  rise  and  the  scan  for  the  pulse  rise  is 
continued.  If  the  auxiliary  condition  is  satisfied  then  point  N  is  considered  to 
be  the  first  point  on  the  pulse  rise.  is  then  set  to  N  -  1. 

Once  has  been  determined  the  next  step  is  to  find  P^,  the  point  number 

for  the  start  of  the  pulse  top.  This  is  done  by  reading  points  P 2*^,  ••• 

from  the  CPUF  and  checking  whether  the  point  P2+I  currently  read  is  significantly 
different  from  the  preceding  point  read.  Thus  SPADE  increments  the  value  of  I 
until  the  condition 

^P^+I  ■  V2+I-I  ^  (6. 2.3.4) 

occurs.  Point  P2+I-I  is  then  possibly  the  first  point  on  the  pulse  top.  To 
check  this  the  values  of  the  next  9  points  must  also  not  differ  from  each  other 
by  more  than  aB.  Thus  if  the  condition 

^P^+I^J  "  ^  (6. 2. 3. 5) 
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holds  for  J  =  1  to  9,  then  is  set  to  P^+I-l. 

Point  P^  is  determined  in  exactly  the  same  way  as  point  P^.  The  top  of  the 
pulse  is  scanned  from  point  P3  and  a  running  average  is  kept.  If 


_P  3+1-1 


(6. 2. 3. 6) 


then  point  P3+I  is  possibly  a  point  0^  the  pulse  drop.  The  next  9  points 


P3+I+I  to  P3+I+9  must  also  differ  from  Yp 


by  at  least  aB/2  according  to  the 


condition  (6.2.3.61.  If  the  next  9  points  do  satisfy  this  condition,  is  set 
to  Pj+I'-l.  P;^  is  then  the  highest  numbered  point  on  the  pulse  top  that  is  not  on 
the  pulse  drop.  However,  if  any  of  the  9  points  do  not  satisfy  this  condition, 
the  search  for  P,  is  continued. 


Point  Pj  is  determined  in  a  similar  way  to  P3.  Thus  SPADE  increments  the 
value  of  I  from  1  until  the  condition 

Yp^+I  -  <  aB  (6.2.3.71 

occurs.  Point  P^+I-1  is  then  possibly  the  first  point  on^  the  background  to  the 
right  of  the  calibration  pulse.  To  check  this,  the  values  of  the  next  9  points 
must  also  not  differ  from  each  other  by  more  than  aB.  Thus  if  the  condition 

^Pi,+I+J  ■  ^P^+I+J-1  ^  (6.2.3.81 

holds  for  J  =  1  to  9  then  P,  is  set  to  P  +1-1. 

5  ‘4 


In  certain  cases,  which  should  rarely  occur,  the  start  or  the  end  of  the 
calibration  pulse  may  be  missing.  SPADE  automatically  checks  for  these 
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conditions.  If  the  start  or  the  end  is  missing,  or  is  set  to  zero, 
respectively,  to  indicate  this  condition. 

SPADE  gives  the  user  the  opportunity  to  re-scan  the  CPDF  several  times  with 
different  values  of  NOISE.  The  first  scan  is  done  automatically,  however,  with 
NOISE  =  10%.  The  pulse  definition  points  found  are  then  printed.  As  illustrated 
in  the  sample  run  in  Section  2.2  of  Reference  5,  several  scans  of  the  data  file 
should  be  made  with  different  values  of  NOISE  (for  example,  20%  and  30%  in 
addition  to  the  one  done  automatically  with  10%).  The  consistency  of  the  results 
should  then  be  checked.  Usually  specifying  higher  values  of  NOISE  is  better  than 
values  less  than  10%. 

When  reasonably  consistent  results  are  obtained,  the  user  should  then  check 
that  the  pulse  definition  points  found  are  correct.  This  can  be  done  either  with 
interactive  graphics  (if  a  graphics  terminal  is  being  used)  or  with  a  special 
listing  routine  in  SPADE  that  allows  a  direct  examination  of  the  data  values  on 
the  CPDF.  In  both  cases  the  pulse  definition  points  can  be  changed  if  required. 
If  interactive  graphics  are  used,  the  pulse  definition  points  can  be  changed  with 
the  graphics  cursor. 

6.2.4  Determination  of  the  Calibration  Pulse  Height 

With  the  pulse  definition  points  defined  on  the  TDF  the  next  step  is  to 
calculate  the  height  of  the  calibration  pulse  in  digitized  units,  Fi(DU),  as 
indicated  in  Figure  5.  This  must  be  done  in  an  accurate  way  that  allows  for 
baseline  drift  and  superimposed  noise  in  each  of  the  regions  defined  by  the  pulse 
definition  points:  Let  P^  to  P^  inclusive  define  Region  1,  P^  to  P,^  inclusive 
Region  2,  and  P^  to  P^  inclusive  Region  3.  Figure  6  shows  these  regions  and 
illustrates  baseline  drift  (although  in  practice  it  would  not  be  as  severe  as 
indicated  and  would  usually  be  undetectable  to  the  naked  eye). 

If  a  least-squares  polynomial  of  low  order  is  fitted  to  the  points  in  the 
regions  1  and  3,  the  result  should  provide  a  reasonable  way  of  smoothing  out  the 
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noise  and  representing  the  unknown  baseline  drift  within  the  calibration  pulse 
itself. 


For  convenience,  the  following  notation  will  be  used.  Let  P  represent  a 
distance  parameter  along  the  CPDF.  If  the  data  points  are  taken  as  separated  by 
one  unit  of  distance  and  the  start  of  the  file  (P  =  0)  is  taken  as  one  distance 
unit  before  the  first  data  point,  then  P  =  1  is  the  distance  of  Point  1  from  the 
start  of  the  file.  P  will  then  be  numerically  equal  to  the  Point  Number  and  will 
measure  the  distance  of  a  data  point  from  the  start  of  the  file.  Non-integral 
values  of  P  will  designate  the  distance  of  positions  between  data  points  from  the 
start  of  the  file.* 

The  least-squares  polynomial  g(P  -  P^)  of  order  m  to  be  determined  using  the 
points  in  Regions  1  and  3  can  be  written 

m+1  .  . 

g(P  -  Pq)  =  C,.(P  -  Pq)^  ^  (0.2. 4.1) 

where  P^  is  a  specific  value  of  P.  P^  is  set  internally  by  SPADE  to  optimise  the 
accuracy  of  the  least-squares  fit  but  can,  in  principle,  be  set  arbitrarily. 

^1"'  ^m+1  coefficients  of  the  polynomial  which  are  determined  from  the 

least  squares  fit  (see  Appendix  C). 

With  g(P  -  P^)  determined,  it  can  be  used  to  estimate  background  value  in 
Region  2.  On  the  pulse  top  in  Region  2,  there  is  a  set  of  discrete  data  points 
that  were  obtained  by  taking  samples  at  equal  time  increments  from  the  pulse  top 
on  the  analog  signal.  Whilst  the  analog  trace  gives  an  accurate  representation 


*  In  the  case  of  a  constant  digitizing  rate,  P  could  be  taken  as  time  measured 
from  the  start  of  a  file.  In  the  case  of  a  variable  digitizing  rate,  however, 
the  time  scale  would  be  nonlinear  with  respect  to  the  P  scale.  It  is  more 
general  and  convenient  to  consider  time  as  a  function  of  P. 
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of  the  signal  for  the  pulse  top,  the  accuracy  of  the  digitized  signal  depends 
upon  the  sampling  frequency  used.  Thus  only  an  approximate  function  can  be 
constructed  to  represent  approximately  the  true  pulse  top.  For  example,  joining 
all  the  digitized  points  along  the  pulse  top  with  straight  lines  would  be  a 
geometric  representation  of  one  possible  approximate  function  for  representing 
the  pulse  top.  Another  would  be  to  fit  a  least-squares  polynomial  to  the  pulse 
top  and  then  to  use  the  polynomial  as  an  approximate  function. 

For  generality,  let  y(P  -  P^)  be  any  approximate  function  for  the  pulse  top. 
Then  h(DU),  the  average  height  of  the  pulse  with  respect  to  the  baseline,  can  be 
determined  from  (see  Figure  6) 


h(DU)  = -  /  [y{P-P„)  -g{P-P„)jdP  (6. 2.4. 2) 

p‘  °  ° 

For  the  evaluation  of  the  integral  in  Equation  6. 2.4. 2,  SPADE  offers  two 
approaches : 

a)  Use  Simpson's  Rule  and  the  known  values  at  the  data  points  on  the  pulse  top. 
The  baseline  data  points  are  evaluated  using  the  polynomial  g(P  -  P^) 
evaluated  at  the  data  points.  The  use  of  this  rule  is  equivalent  to 
constructing  the  approximate  function  y(P  -  P^)  with  quadratic  polynomials 
as  interpolation  functions  over  consecutive  groups  of  three  data  points. 

This  method  will  be  clarified  in  Section  6.7  where  Simpson's  Rule  is 

di scussed. 

b)  Fit  a  least-squares  polynomial  to  the  points  on  the  pulse  top  and  then  use 
this  polynomial  as  the  approximate  function  y(P  -  P^).  Equation  (6. 2.4. 2) 
can  then  be  integrated  in  closed  form  to  find  the  average  pulse  height 
h(DU). 

Choice  a)  will  usually  be  the  most  accurate  and  should  be  chosen  in  most 
cases.  Specifying  various  orders  of  polynomial  for  the  background  will  then  give 
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an  estimate  of  accuracy.  Choice  b)  can  be  selected  as  a  check  if  required.  It 
is  important  to  note,  however,  that  the  last  value  of  h(DU)  calculated  is  the  one 
that  will  be  used.  Thus,  if  checks  on  h(DU)  are  carried  out,  the  one  required 
should  be  re-calculated  before  leaving  the  calculation.  Once  a  satisfactory 
value  of  h(DU)  has  been  calculated,  the  Scale  Factor  E  is  calculated  in 
accordance  with  Equation  (6.2. 1.1). 

6.3  Option  3:  Scale  the  Transducer  Signal  Data 

6.3.1  Introduction 


The  object  of  Option  3  is  to  convert  the  data  on  the  SDF  for  a  specified 

transducer  from  digitized  units  to  real  units  by  applying  the  value  of  the  Scale 

Factor  E  computed  with  Option  2.  This  value  of  E  is  stored  on  the  TDF.  SPADE 
will  check  that  E  for  the  specified  transducer  has,  in  fact,  been  calculated.  If 
it  has  not,  then  an  information  message  will  be  given  and  Option  2  will  be 
required  to  be  done  before  Option  3  can  be  executed. 

At  the  start  of  Option  3,  SPADE  will  list  the  transducers  that  have  been 
defined  on  the  TDF.  The  Channel  and  Transducer  Code  will  then  be  requested  to 
specify  the  SDF  that  is  to  be  scaled.  A  full  illustration  of  the  use  of  this 

option  is  given  in  the  sample  run  in  Section  2.3  of  Reference  5. 

6.3.2  Definition  of  the  Signal  Definition  Points 

Figure  7  shows  a  schematic  illustration  of  a  digitized  transducer  signal 
with  superimposed  background  drift  and  noise.  Typically  the  total  number  of 
points  on  the  SDF  is  5000.  As  is  done  for  the  CPDF,  SPADE  uses  a  set  of 
transducer  Signal  Definition  Points  (SDP)  denoted  as  P,,  P.,  P,,  P  ,  P,,  and  P,. 
These  points  are  stored  on  the  TDF.  P^  to  P^  will  now  be  defined. 

Points  Pj  to  P^  inclusive  will  be  defined  by  SPADE  as  a  set  of  points  before 
the  signal  that  is  used  for  background  correction  (Region  1).  Figure  7 
illustrates  this  definition.  Usually  Pj  will  be  set  to  1  to  denote  the  first 
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point  on  the  SDF.  Also  will  be  set  to  the  point  number  that  defines  the  start 
of  the  signal . 

Points  Pj  and  Pj^  are  used  only  for  pressure  transducer  signals.  They  are 
not  used  for  other  types  of  transducer  signal  in  the  current  version  of  SPADE. 

For  a  pressure  signal  P3  is  set  to  the  Point  Number  of  the  maximum  pressure  and 
Pi^  to  the  Point  Number  that  defines  the  end  of  the  positive  phase.  Figure  8 
shows  a  typical  type  of  pressure  pulse  with  P3  and  P^  defined. 

Points  Pj  to  Pg  inclusive  define  a  set  of  points  after  the  signal  (Region  3) 
that  is  used  for  background  correction.  Usually  P^  is  be  set  to  0  to  indicate 
to  SPADE  that  no  points  after  the  signal  are  to  be  used  for  background 
correction. 

Points  P^  to  P^  exclusive  will  define  Region  2.  This  region  contains  the 
actual  signal.  If  P^  is  zero,  Region  3  does  not  exist  and  the  end  of  Region  2  is 
undefined.  For  acceleration  and  strain  this  is  generally  the  case  since  some 
signal  is  often  detectable  at  the  end  of  the  SDF.  For  pressure  signals,  it  may 
be  possible  to  define  P^. 

6.3.3  Determination  of  the  Signal  Definition  Points 

P^  is  set  to  1  and  will  remain  at  this  value  unless  explicity  changed  by  the 
user.  SPADE  will  automatically  scan  the  SDF  for  the  start  of  the  signal  and  set 
P^  to  the  PN  found.  The  technique  used  to  scan  for  P^  is  precisely  the  same  as 
that  used  to  scan  for  P^  in  the  case  of  the  Calibration  Pulse.  This  technique 
was  discussed  in  Section  6.2.3.  A  parameter  NOISE  is  again  used.  The  value  of 
NOISE  is  initially  set  to  10/^  but  several  values  should  be  tried  until  consistent 
values  for  P^  are  obtained.  For  pressure  signals  which  have  steep  rise  times 
(for  pressure  transducers  on  the  top  surface),  values  of  say  5«  and  15%  may  give 
results  which  are  consistent  with  10%.  For  strain  gauges  and  accelerometers,  it 
may  turn  out  that  values  from  1%  to  5%  will  give  consistent  results.  The  lower 
noise  level  factor  that  may  be  required  in  the  latter  case  for  consistent  results 
is  due  to  a  slower  rise  time. 
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When  P2  has  been  determined,  the  user  is  given  a  choice  of  checking  the 
results  using  interactive  graphics  or  a  special  listing  routine  that  allows 
direct  examination  of  the  data  values  on  the  SDF.  At  this  stage  SPADE  also 
allows  the  user  to  redefine  the  definition  points  Pj  and  P^.  Points  P3  and  P^ 
only  apply  to  pressure  pulses  and  are  determined  automatically  with  Option  5. 

Thus  they  need  not  be  manually  set  within  this  option  although  they  could  be  if 
desired. 

Points  P^  and  P^  will  usually  not  be  required.  The  condition  P^  =  0,  which 
is  the  default  value,  will  flag  the  condition  for  no  background  correction  after 
the  signal.  However,  SPADE  will  ask  whether  background  correction  is  to  be 

applied  using  points  after  the  signal.  If  "Yes"  is  indicated,  P5  and  P^  will 
have  to  be  given  for  defining  Region  3. 

6.3.4  Background  Correction  and  Scaling  of  the  Signal  Data  File 

With  the  SOP  defined  for  the  SDF,  the  next  step  is  to  correct  for  background 
and  scale  the  data  to  real  units.  The  processed  data  will  then  be  placed  in  the 
Scaled  Signal  Data  File  (SSDF). 

Figure  7  shows  a  schematic  illustration  of  a  transducer  signal  with  the 
background  drift  and  noise  exaggerated.  To  correct  for  background,  the  same 
procedure  as  described  in  Section  6.2,4  for  the  CP  is  followed.  For  the  case  of 
a  transducer  signal,  however,  only  the  points  in  Region  1  (Figure  7)  are  usually 
used  for  background  correction.  The  least-squares  polynomial  is  then 
extrapolated  into  the  signal.  Let  g(P  -  P^)  be  the  polynomial  as  defined  in 
Equation  (6. 2.4.1).  Also  let  Yp  be  the  digitized  data  value  of  Point  P.  Then 
the  data  values  in  real  units,  denoted  by  Zp,  are  given  by 

Zp  =  E|Vp  -  giP-Pp)  ]  (6.3.4. 1) 

where  E  is  the  scale  factor  determined  in  Option  1  or  Option  2. 

The  real  data  values,  Zp,  are  written  to  an  SSDF  created  by  SPADE  with  a 
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natne  given  according  to  the  file  naming  convention  discussed  in  Appendix  B.  The 
SDF  is  not  overwritten. 


6.4  Option  4:  Determine  the  Start  Time  for  the  Pressure  Transducer  Signals 


Each  point  on  the  CPDF  corresponds  to  a  definite  time  from  the  IRIG  B  signal 
which  provides  a  common  time  base  for  all  events.  Let  measure  distance  along 
the  SDF  or  the  SSDF  for  the  i'th  transducer  on  the  TDF  (recall  that,  by  defini¬ 
tion,  p’  will  be  the  actual  Point  Number  at  a  data  point  since  all  data  points 
are  considered  to  be  equally  spaced  by  1  unit  of  distance  and  p’=  0  is  1  distance 
unit  after  the  first  point).  Also  let  T^  be  the  start  time  of  digitizing  (from 
the  IRIG  B  signal)  for  transducer  i.  With  the  digitizing  rate  R  ,  the  time  T 
corresponding  to  P’  can  be  found  from  the  following  relation  (see  Figure  9) 


T  =  t’  +  pVr^ 


(6.4.1) 


In  general,  it  is  not  convenient  to  use  the  IRIG  B  time  T  but  rather  a  new 
time  scale  t  defined  as 


^  ^  ■  ^0  (6.4.2) 

where  T^  is  an  IRIG  3  time  chosen  as  the  t  =  0  point  on  the  new  time  scale.  No 
superscript  i  is  required  in  Equation  (6.4.2)  since  the  IRIG  B  time  scale  and  the 
new  time  scale  are  common  time  scales  for  all  transducers.  T^  will  correspond  to 
P^  from  the  file  start  for  the  i'th  transducer,  which  will  usually,  although  not 
necessarily,  fall  within  the  file  limits  on  Point  Number.  For  transducer  i,  for 
example,  we  have  from  Equation  (6.4.1) 

To  =  t’  4  P’/R^  (6.4.3) 

Thus  once  a  T  has  been  selected  on  the  IRIG  B  time  scale,  P^  can  be  calculated 
0  .  ’0 

by  solving  6.4.3  for  P^  as 
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(6.4.4) 


It  is  important  to  note  that,  in  general,  t’  will  be  different  for  each 
transducer  i. 

can  essentially  be  chosen  arbitrarily;  the  choice  simply  amounts  to 
choosing  a  common  time  point  for  t  =  0  for  all  transducers.  SPADE  gives  two 
choices  for  choosing  T^.  These  choices  are  described  as  choice  I  and  II  below. 


I.  Choose  Tq  to  be  the  average  IRI6  B  time  of  arrival  of  the  pressure 
pulses  from  all  the  pressure  transducers  on  the  top  surface  of  the 
panel.  This  choice  is  appropriate  for  a  symmetrically  placed  charge  as 
shown  in  Figure  2. 

II.  Choose  Tq  as  the  earliest  IRIG  8  time  of  arrival  of  the  pressure  pulses 
from  a  pressure  transducer  on  the  top  panel  surface.  This  choice  would 
be  appropriate  if  the  charge  over  the  plate  in  Figure  2  was  not 
symmetrically  placed. 


The  choice  of  I  or  II  should  be  made  for  convenience  only.  Pressure 
transducers  on  the  top  surface  are  convenient  to  use  in  both  cases  since  they  are 
directly  exposed  to  a  blast  wave  and  thus  have  very  steep  rise  times. 

The  procedure  followed  by  SPADE  will  now  be  described.  The  overall  purpose 
of  Option  4  is  to  determine  a  P^  for  each  transducer  i  on  a  TDF  corresponding  to 
a  specific  Series,  Shot,  and  Panel  and  to  write  the  values  to  the  TDF.  These 
values  can  then  be  used  by  other  options  of  SPADE  and  by  SPADE-PLOT  |4  to 
establish  an  appropriate  time  base.  For  the  purpose  of  the  following  discussion, 
it  will  be  assumed  that  transducers  i  =  1  to  Np  are  pressure  transducers  on  the 
top  panel  surface.  The  total  number  of  transducers  on  the  TDF  is  N. 

1.  Extract  the  start  of  digitizing  time  t’  for  i  =  1  to  N  from  the  TDF. 

2.  Determine  the  minimum  digitizing  start  time  T^j  for  the  pressure  trans- 
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ducers  on  the  top  surface.  Thus  =  minimum  of  t’  for  i  =  1  to  Np. 

3.  For  all  transducers  i  =  i  to  N 

Set  t’  — ^  (6.4.5) 

S  S  M 

This  just  shifts  the  zero  point  on  the  IRIG  8  time  scale  for  computational 
convenience.  In  steps  4  to  6,  reference  to  IRIG  B  time  will  be  assumed 
measured  from  this  new  starting  point. 

4.  Calculate,  for  i  =  1  to  Np  from  Equation  (6.4.1) 

'  ♦ 

t]  =  t1  +  p’/r’  (6.4.6) 

where  P2  is  the  second  SOP  for  the  i'th  pressure  transducer,  t’  is  then  the 
IRIG  B  time  (measured  from  the  minimum  digitizing  start  time)  for  the  start 
of  the  pressure  pulse. 


5.  If  choice  I  is  being  considered  for  determining  T^,  set 


T 


0 


P 


(6.4.7) 


Tq  is  then  the  average  time  of  arrival  of  the  pressure  pulses. 

If  choice  II  is  being  considered,  set  T^  to  the  minimum  of  the  T^ 
values.  T^  will  then  be  the  first  time  of  arrival  of  a  pressure  pulse. 

6.  For  i  =  1  to  N  calculate  the  P^  values  for  all  transducers  from 
Equation  (6.4.4).  The  values  of  P^  then  define  the  Point  Numbers  on  the  SDF 
and  SSDF  that  correspond  to  the  t  =  0  point. 
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6.5  Option  5:  Determine  Pressure  Signal  Characteristics 
6.5.1  Introduction 


Figure  10  shows  a  typical  pressure  pulse  as  would  be  detected  by  a  pressure 
gauge  exposed  directly  to  a  blast  wave.  The  curve  represents  a  plot  of 
overpressure  against  time  t.  It  is  the  purpose  of  this  option  to  determine  the 
pressure  pulse  characteristics  that  are  indicated  in  Figure  10. 

Section  6.5.2  will  first  review  the  definitions  of  the  pressure  pulse 
characteristics  indicated  in  Figure  10.  Section  6.5.3  will  present  various 
functional  relations  that  can  be  used  to  approximate  pressure  pulses.  Finally, 
Sections  6.5.4  and  6.5.5  will  discuss  how  SPADE  uses  the  various  functional 
relations  to  describe  pressure  signals. 

6.5.2  Definition  of  Pressure  Pulse  Characteristics 

The  pressure  pulse  illustrated  schematically  in  Figure  10  is  of  typical 

shape  (Reference  6,  page  111  or  Reference  7,  pages  4-8)  and  represents  an 

undisturbed  or  side-on  shock  wave  as  it  propagates  through  the  air.  The  pressure 

pulse,  when  detected  by  a  pressure  guage,  will  have  a  certain  time  of  arrival  T  . 

d 

For  T  <  T,  the  pressure  is  zero  relative  to  ambient  pressure.  For  T  >  T  a  steep 

a  d 

rise  in  pressure  occurs  to  a  maximum  at  time  T^  +  T^  known  as  the  peak 

overpressure  (or  peak  side-on  overpressure)  P^.  T^  is  known  as  the  rise  time 

and  is  assumed  to  be  zero  for  an  ideal  blast  wave.  After  the  maximum  has  been 

reached,  the  pressure  pulse  decays  to  zero  and  then  becomes  negative  to  a  maximum 

negative  amplitude  of  P  known  as  the  peak  underpressure.  If  the  time  at  which 

^  +  + 

the  pressure  pulse  becomes  negative  is  T^+T  ,  then  T  is  referred  to  as  the 

+  ^  +  - 

positive  phase  duration.  For  T  +T  «;  T  <  T  +T  +T  the  negative  phase  (or  suction 
phase)  occurs;  that  is,  a  partial  vacuum  exists  behind  the  shock  front.  T”  is 
known  as  the  negative  phase  duration.  For  T  .■  T^+T^+T"  the  pressure  is  assumed 
to  remain  at  zero. 
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The  positive  specific  impulse,  i^,  is  the  area  under  the  positive  phase  of 
the  curve  and  represents  the  impulse  per  unit  area.  Thus,  i^  is  defined  as 

ij  =  j  P(t)dt  (6. 5. 2.1) 

a 

and  i^,  the  negative  specific  impulse,  is  defined  as 

Ta+T'^+T" 

i"  =  -/  P(t)dt  (6. 5. 2. 2) 

T  +T‘^ 
a 

SPADE  calculates  the  following  pressure  signal  characteristics  directly  from 

tfie  experimental  data:  the  time  of  arrival,  peak  overpressure,  rise  time, 

positive  phase  duration  and  positive  specific  impulse.  In  addition,  the  signal 

definition  points  P3  and  P^^  are  calculated  to  indicate  the  time  of  peak 

overpressure  and  the  end  of  the  positive  phase,  respectively.  These  two 

additional  points  are  written  to  the  TOF  and  allow  T  and  T^  to  be  calculated 

r 

when  required.  In  addition,  the  pressure  SSDF  is  integrated  to  give  a 
corresponding  file  containing  specific  impulse  against  time  from  which  the 
maximum  specific  impulse  and  the  time  of  occurrence  can  be  determined. 

6.5.3  Functional  Forms  Used  to  Approximate  Pressure  Signals 

The  following  functional  forms  for  representing  pressure  pulses  are  taken 
from  Reference  7. 

For  the  positive  phase  only 

P(t')  =  P^  exp  (-b^t'/T^)  (6. 5. 3.1) 
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where  is  a  constant  known  as  the  decay  constant,  t'  is  the  time  measured  from 
T^+T^  (thus  t'=T-T^-T^).  This  function  can  be  used  to  obtain  a  reasonably  good 
fit  usually  over  at  least  50%  of  the  positive  phase. 

For  both  the  positive  and  negative  phase  a  very  crude  linear  approximation 
i  s 


P(t')  =  pJ(l-t7T^)  (6. 5. 3. 2) 

If  Equation  (6. 5. 3. 2)  is  combined  with  Equation  (6.5.3. 1)  however,  a  much 
improved  relation  results  as  follows: 

P(t')  =  P^d-t'/T"^)  exp  (-b^t'/T'^)  (6. 5. 3. 3) 

Equation  (6. 5. 3. 3)  is  sometimes  referred  to  as  the  modified  Friedlander 
equation  [6].  This  functional  form  will  usually  give  a  reasonable  fit  to 
experimental  data  up  to  a  time  T=T,+3T^/2. 

a 

In  Equation  (6. 5. 3. 3)  the  decay  constant  is  assumed  to  be  constant.  An 
improvement  is  obtained  if  a  time  dependent  decay  constant  b  is  introduced  and 
assumed  to  be  a  function  of  time  t'  as 


b 


1  + 


l+hli,:) 


(6. 5. 3. 4) 


where  g  and  h  are  constants.  A  somewhat  simpler  but  more  restrictive  form  is 
sometimes  adequate  as  follows: 


b  -  b  [1  -  f(^)] 
°  T 


(6. 5. 3. 5) 
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where  f  is  a  constant. 

An  equation  for  P(t')  that  is  known  to  give  an  excellent  fit  to  experimental 
blast  waves  is  [7] 


**  "  /T^  "  /T^ 

P(t')  =  P^[l  -(^)J|ae~“^  +  (l-a)e‘^^  \  (6. 5. 3. 6) 

where  a,  a  and  3  are  constants. 

Whilst  Equations  {6.5.3.31  and  (6. 5. 3. 6)  are  valid  for  both  the  positive  and 
negative  phases,  the  negative  phase  may  not  be  fitted  with  good  accuracy.  An 
equation  that  is  known  to  give  a  good  fit  to  the  negative  phase  is 

P(t)  =  -p"(^^)(l  -  iiL)exp[-4^-^h  (6.5.3.71 

T"  T"  T' 

where  t  is  now  measured  from  T^. 

6.5.4  Direct  Calculation  of  Pressure  Pulse  Parameters 


All  of  the  equations  (6.5.3.11  to  (6.5.3.31  and  (6.5.3.61  involve  P^  and  T^. 
Any  remaining  constants  can  be  determined  by  using  the  experimental  data  and  the 
appropriate  equations.  For  example,  Equation  (6.5.3.21  can  be  integrated  to  give 
an  expression  for  i^.  in  the  equation  can  then  be  determined  using  the  i^ 
found  by  integration  of  the  experimental  values.  Thus  if  P(t'l  given  by 
(6.5.3.21  is  substituted  into  Equation  (6.5.2.11  and  the  resulting  equation  for 
i^  is  solved  for  there  results 


T 


+ 


(6.5.4.11 
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The  experimental  value  of  found  by  integration  of  the  pressure  time  data 
can  then  be  substituted  into  Equation  (6. 5.4.1)  together  with  an  experimental 
value  of  to  give  T^. 

A  similar  procedure  can  be  followed  for  Equation  (f.5.3.1).  If  i^  is 
calculated  using  Equation  (6.5.3. 1),  the  following  equation  is  obtained. 


exp(-b^)  + 


.+ 


s 


(6. 5.4. 2) 


In  this  case  i^,  and  could  be  supplied  and  the  nonlinear 
Equation  (6.5.4. 2)  solved  iteratively  for  b^.  This  value  of  b^  would  then 
guarantee  that  i^  calculated  with  Equation  (6.5.3. 1)  be  the  same  as  the 
experimental  value. 


In  the  case  of  Equation  (6. 5. 3. 3)  a  more  complicated  relation  results. 

P  T  -  1b 

exp(-b^)  +  ■  -V  T-  ^  1  (6. 5. 4. 3) 

U  pTjT  y 

S 

Again,  i^,  P^  and  T^  could  be  supplied  to  calculate  a  b^.  Other 
possibilities  would  be  to  also  calculate  i^  and  then  to  determine  T^  to  match  the 
experimental  i~. 

SPADE  offers  facilities  for  calculating  the  constants  in  equations  (6. 5.3.1) 
to  (6. 5. 3. 3)  using  the  procedures  mentioned.  T^  can  be  determined  from 
(6.5.4. 1).  i^  and  P^  are  calculated  from  the  experimental  data  on  the  SSDF.  b^ 
is  calculated  from  equations  (6. 5.4.2)  or  (6. 5. 4. 3)  using  experimental  values  of 
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6.5.5  Least  Squares  Calculation  of  Parameters 


SPADE  also  allows  the  parameters  in  Equations  (6. 5. 3.1)  and  (6. 5. 3. 3)  to  be 
calculated  using  least  squares.  The  advantage  of  this  technique  is  that  all 
experimental  data  points  can  be  used  to  estimate  the  parameters.  The  approach 
taken  will  now  be  outlined  for  the  case  of  Equation  (6. 5. 3.1). 


Equation  (6.5.3. 1)  is  repeated  below  for  convenience. 


P(t')  =  P^exp(-bQ^) 


Taking  the  natural  logarithm  of  both  sides  then  gives 


ln(P)  =  In(P^)  -  b^^ 


(6.5.3. 1) 


(6. 5.5.1) 


Values  of  P  for  each  time  point  on  the  SSDF  are  available.  Thus  if  a  linear 

least  squares  fit  for  points  Ct',  In  (P)]  is  done  over  a  selected  region  of  the 

positive  duration,  the  slope  of  the  straight  line  obtained  will  give  b  /T^  and 

the  ordinate  intercept  will  give  ln(P  ).  Thus,  provided  a  T  is  available  or  can 

^  + 

be  estimated,  b^  can  be  found  and,  of  course,  P^  can  be  found  by  taking  the 
exponential  of  In(P^). 

If  the  parameters  are  determined  by  this  procedure,  SPADE  will  request 
values  of  and  E^  to  define  the  fraction  of  the  time  domain  (Ta+Tr  to  Ta+T^)  to 
be  used  to  do  the  least  squares  fit.  Only  points  on  the  pressure  SSDF  with  PN's 
between  E,(P  -PJ  +  P,  and  E  ,(P  -P  )  +  P,,  inclusive,  will  be  used  in  the  least 
squares  fit.  As  indicated  in  Section  6.5.2,  the  SDP  P^  and  P^  identify  the 
position  of  maximum  pressure  and  positive  duration,  respectively. 
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The  calculation  for  and  by  this  procedure  is  outlined  in  the  sample 
run  in  Section  2.5  of  Reference  5.  Interactive  graphics  are  very  useful  in  this 
option  for  determining  whether  the  values  found  are  satisfactory.  Figure  11 
illustrates  a  series  of  plots  obtained  with  interactive  graphics.  Figures  11(a) 
to  11(c)  indicate  a  fit  to  a  pressure  pulse  with  =0,  and  =  0.3  which  are 
typical  and  the  suggested  (default)  values.  There  are  about  100  points  in  the 
positive  phase.  An  excellent  fit  is  obtained  over  about  one  half  the  positive 
phase.  The  result  over  the  whole  positive  phase  is  quite  good  in  fact. 

The  interactive  graphics  procedure  allows  "windowing"  to  be  done  with  the 
graphics  cursor.  The  small  box  indicates  the  portion  to  be  windowed. 

Figure  11(c)  indicates  the  result  after  windowing  to  the  limit  of  digitizing. 

The  digitized  points  can  be  clearly  seen  joined  with  straight  lines. 

Figure  11(d)  shows  the  result  of  using  E^  =  0,  E_,  =  0.8.  Because  points  in 
the  latter  part  of  the  positive  duration  are  used,  a  loss  of  accuracy  occurs  for 
point  numbers  0  to  20.  This  result  would  thus  not  be  acceptable  since  a  good  fit 
should  be  obtained  for  at  least  the  first  one  third  of  the  positive  duration. 

6.6  Option  6;  Integrate  Acceleration  to  Obtain  Relative  Velocity 

6.6.1  Introduction 


This  option  requires  very  little  interaction  of  the  user  with  SPADE.  An 
SSDF  containing  accelerometer  output  can  be  integrated  in  time  to  obtain  an  SSOF 
containing  velocity.  Simpson's  integration  technique  [8]  is  used  but  not  in  the 
standard  way.  The  required  formulas  will  be  developed  in  this  section. 

6.6.2  Theory 

To  find  velocity  it  is  necessary  to  integrate  acceleration  with  respect  to 
time  t.  If  R  is  the  digitizing  rate,  there  are  1/R  (=h)  time  units  between  each 
sample  on  the  SSOF.  Let  be  the  point  number  on  the  SSDF  corresponding  to 
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t  =  0.  Then  the  time  t  at  point  P  is 

t  .  h(P-P„)  (6.6.2.1) 

Let  y(t)  represent  the  true  acceleration  at  a  time  t  on  the  SSDF  containing 
accelerations.  Because  of  the  finite  sampling  rate  in  the  digitizing  process 
y(t)  is  only  known  when  t  =  t^.,  where  t^.  is  the  time  corresponding  to  the  i'th 
point  number.  To  construct  an  approximate  acceleration  function,  quadratic 
interpolation  will  used  to  give  values  of  acceleration  between  sample  points. 
Figure  12(a)  illustrates  this  approximation  procedure.  The  SSDF  is  assumed  to 
contain  only  seven  points  in  this  case.  The  solid  curve  represents  the  true  y(t) 
which  is  only  actually  known  at  the  data  points.  The  dotted  curve  represents 
y(t),  an  approximate  function.  Points  1  to  3  have  a  single  quadratic  polynomial 
through  them  (3  points  will  determine  a  quadratic  polynomial  uniquely).  Points  3 
to  5  have  a  second  quadratic  througli  them  and  Points  5  to  7  a  third.  The  ♦'irst 
and  second  quadratic  have  a  common  value  at  Point  3  but  the  slopes  do  not  match. 
The  same  situation  occurs  at  Point  5.  y(t)  is  thus  a  continuous  function  but  is 
only  piecewise  differentiable  since  the  slope  is  discontinuous  at  points  where 
the  quadratics  are  joined.  y(t)  can  thus  be  integrated  in  a  piecewise  fashion  to 
provide  an  approximation  to  the  integation  of  the  true  acceleration  y(t). 

To  develop  the  formulas  necessary  for  integration,  consider  the  region  in 
Figure  12(b)  defined  by  the  points  corresponding  to  times  t^,  t^+h,  and  tj+2h. 

Let  the  quadratic  polynomial  defining  the  acceleration  over  this  region  be 
defined  as 


yj(t)  =  b^+b2(t-t^)+b3(t-tl)^  (6.6. 2. 2) 

The  coefficients  b^,  b^  and  b^  in  Equation  (6. 6. 2. 2)  can  be  found  by  substituting 
t  =  tj,  t  =  t^  +  h  (or  t^),  and  t  =  tj  +  2h  (or  tj)  into  Equation  (6. 6. 2. 2).  The 
conditions  yi(t^)  =  yj,  yi(t2)  =  y^,  and  yi(tj)  =  y where  y^  y^,  and  yj  are 
known  values  of  accleration,  determine  the  coefficients  b^,  b^,  and  b^  uniquely 

I 
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(6.6. 2.3) 

_  4y2  -  3yi  -  yj 

2"  (6. 6. 2. 4) 

yi  -  2y2  +  yj 

b3  =  -  (6.6. 2.5) 

2h^ 


The  next  step  is  to  derive  an  expression  for  yj(t)  integrated  from  t^  to 
t^+h  and  t^+2h.  For  generality  let 


tj  +  ah 
^1 


tj+oh 

=  /  y(t)dt 


(6.6. 2.6) 


where  a  is  a  fraction  of  the  time  h  between  two  data  points.  If  Equation 

(6. 6. 2. 2)  is  substituted  into  (6.6. 2.6)  and  the  equations  (6. 6. 2. 3)  to  (6.6. 2. 5) 

are  utilized,  one  can  show  after  a  straightforward  reduction  that 


ti+ah  312  2 

I,.  =  had - a  +  —  a  )yi  "*■  ha  (1 

4  6 

Thus,  setting  a  =  1  gives 


1  ^  u  lx 

-  aly,  +  ha  (-  a  —  )y  3 

3  6  4 


(6.6. 2.7) 


tj+h 

d  =  — (5yi  +  8y2  -  y3)  (6.6. 2.8) 

d  12 


and  setting  a  =  2  gives 


ti+2h  , 

I  =  -  (yi  +  4y2  +  ya)  (6.6. 2.9) 

3 
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Thus 


ti+2h  t^+2h  tj+h  . 

+h  =  =  —  (-yi+  8y2  +  Syj)  (6.6.2.10) 


Equations  (6. 6. 2. 8)  and  (6.6.2.10)  are  the  equations  utilized  by  SPADE  to 
integrate  acceleration  to  obtain  velocity.  Equation  (6.6. 2.8)  allows  an 
integration  over  the  first  half  of  the  interval  shown  in  Figure  12(b)  and 
(6.6.2.10)  over  the  second  half.  Thus  by  integrating  over  each  individual 
segment  in  each  of  the  three  regions  in  Figure  12(a)  the  velocity  at  each  time 
point  can  be  found.  Thus 


ti+h 

^2  =  If  =  —  (5yi  +  8y2  -  y^) 

i2 

t^+2h  tj+h  . 

V3  =  I.  =1.  +  —  (-yi  +  8y2  +  Sy^) 

‘■i  ‘■l  12 

t,+3h  tj+2h  . 

V.  =  L  =4  +  ^  (Syj  +  8y,  -  y^) 

12 


ti+4h  t+3h  . 

V5  =  4  =  It  +  —  (-ys  +  8y4  +  Sy^) 

12 


Thus  velocities  at  all  the  discrete  time  points  can  be  obtained. 

It  is  perhaps  worth  mentioning  that  Equation  (6. 6. 2. 9)  is  usually  the  way 
that  the  Simpson's  Rule  is  presented;  that  is,  a  in  Equation  (6.6. 2.7)  is 
restricted  to  being  an  even  number.  The  advantage  of  using  the  more  general  form 
in  Equation  (6.6. 2.7)  is  that  a  can  be  chosen  to  be  any  value  from  0  to  2 
inclusive.  Thus,  if  required,  any  number  of  points  within  each  region  can  be 
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chosen  at  which  the  velocity  is  required.  As  indicated  above,  however,  SPADE 
only  evaluates  the  velocities  at  the  same  time  points  for  which  the  acceleration 
is  defined. 

6.6.3  SPADE  Processing  for  Velocity 

SPADE  will  initially  request  the  Series,  Panel,  and  Shot  numbers  as  usual. 
After  Option  6  has  been  flagged,  SPADE  will  request  the  Channel  Number  and 
Transducer  Code.  If  the  Transducer  Code  does  not  correspond  to  an  accelerometer, 
SPADE  will  request  the  Channel  Number  and  Transducer  Code  to  be  respecified. 

Also,  since  the  integation  is  done  from  the  t  =  0  position  on  the  SSDF  (time  T^ 
on  IRIG  B),  SPADE  will  check  that  P^  has  been  determined  with  Option  4  (see 
Section  6.4).  On  completion  of  the  integration,  SPADE  will  indicate  that  the 
file  containing  the  velocities  has  been  created.  The  original  acceleration  SSDF 
is  not  overwritten;  a  new  file  referred  to  as  a  Derived  Velocity  Signal  File 
(DVSF)  is  created  with  a  file  name  according  to  the  convention  in  Appendix  B2. 

The  format  of  the  DVSF  is  identical  to  that  of  the  corresponding  SSDF  containing 
the  accelerations.  Thus  if  on  the  SSDF  is  the  acceleration  at  the  N'th  PN, 
then  on  the  DVSF  it  is  the  velocity  at  the  same  PN. 

6.7  Option  7:  Integrate  Velocity  to  Obtain  Relative  Displacement 

This  option  is  esentially  identical  to  Option  6.  In  this  case,  however, 
there  are  two  choices:  1)  To  integrate  an  SSDF  containing  data  from  a  velocity 

transducer  or  2)  To  integrate  a  DVSF.  In  each  case,  SPADE  will  check  that  Option 

4  has  been  completed  so  that  integration  can  be  carried  out  from  the  t  =  0 
position.  On  completion  of  the  integration,  SPADE  will  indicate  that  the  file 
containing  the  relative  displacements  has  been  created.  The  original  velocity 
SSDF  (choice  1)  or  DVSF  (choice  2)  is  not  overwritten;  a  new  file  referred  to  as 
a  Derived  Displacement  Signal  File  (DDSF)  is  created.  The  file  naming  convention 

discussed  in  Appendix  B2  will  distinguish  a  DDSF  as  having  been  created  from  an 

SSDF  containing  velocities  or  a  DVSF.  The  format  of  the  DDSF  is  identical  to 
that  of  the  corresponding  SSDF  or  DVSF.  Thus  if  on  the  SSDF  or  DVSF  is  the 
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velocity  at  the  N'th  point  number,  then  on  the  DDSF  it  is  the  displacement  at  the 
same  time  point.  It  is  worth  noting  that  displacement  data  may  be  recorded 
during  the  tests  using  displacement  measuring  transducers.  This  displacement 
data  would  be  placed  on  a  SSDF  according  to  the  above  file  format  and  the  file 
would  be  named  by  the  convention  presented  in  Appendix  B2. 

6.8  Option  8:  Obtain  Acceleration,  Velocity  and  Displacement 
Based  on  Fit  to  Displacement 

6.8.1  Introduction 


In  the  execution  of  Option  6  and  Option  7,  several  effects  are  produced  due 
to  background  noise  effects  and  possibly  unknown  initial  conditions.  For 
example,  if  a  signal  constant  in  time  is  superimposed  on  an  acceleration  signal, 
the  velocity  signal  will  have  a  signal  linear  in  time  superimposed  on  it. 
Similarly,  because  of  double  integration,  the  displacement  signal  will  have  a 
signal  quadratic  in  time  superimposed  on  it.  Although  Option  3  will  correct  for 
background  to  a  significant  extent,  small  non-zero  acceleration  baseline  effects 
may  still  be  present  and  turn  out  to  be  significant. 

Another  effect  that  can  cause  extraneous  results  is  unknown  initial 
conditions.  An  acceleration  signal  from  an  accelerometer  can  only  be  used  to 
extract  information  on  changes  of  velocity.  After  execution  of  Option  6,  the 
resulting  DVSF  will  only  contain  relative  velocities.  Unless  the  initial 
velocity  is  zero,  the  Derived  Velocity  Signal  (DVS)  will  require  a  signal 
constant  in  time  to  be  superimposed  on  it  to  give  the  correct  initial  velocity. 
This  will  then  cause  a  signal  linear  in  time  to  be  superimposed  on  the  Derived 
Displacement  Signal  (DD5). 

In  Section  6.8.2,  the  above  considerations  will  be  presented  in  more  detail 
together  with  the  method  used  by  SPADE  to  obtain  corrected  accelerations, 
velocity,  and  displacement  signals.  Section  6.8.3  will  then  describe  how  SPADE 
proceeds  with  Option  8. 
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6.8.2  Theory  and  Method  of  Correction 

Before  proceeding  with  a  description  of  this  option,  it  will  be  useful  to 
define  an  integration  operator,  I^.  If  f(t)  is  a  function  of  time,  ij  is  defined 
as 

ijflt)  =  /Jflxldx  (6.8. 2.1) 

where  x  is  a  dummy  variable  of  integration.  The  operation  can  also  be  repeated 
as  follows: 


ijljf(t)  =  /J  /^f(x)dxdy  (6.8. 2.2) 

y  is  another  dummy  variable  of  integration. 

Now  let  a(t)  be  the  true  acceleration  time  signal  for  a  point  on  a 
structure.  Also  let  a(t)  be  the  measured  value  of  acceleration.  Then  the 
measured  acceleration  time  signal  a(t)  can  be  regarded  as  the  true  acceleration 
time  signal  with  a  superimposed  error  signal  c  (t).  Thus 

a(t)  =  a(t)  +  'fit)  (6.8. 2.3) 

When  Option  6  is  carried  out,  the  SSDF  containing  sampled  values  of  a(t)  is  used 
to  create  a  DVSF  containing  values  of  v(t),  the  derived  relative  velocities. 

From  Equation  (6.8. 2. 3)  we  can  write 

v(t)  =  Iq  a(t)  =  ija(t)  +  ij  e^t)  (6.8. 2.4) 

Now  I^a(t)  =  v(t)  -  v(o)  where  v(t)  is  the  true  velocity  obtained  from  a(t)  by 
integration  and  v(o)  is  the  true  value  of  velocity  at  t  =  0.  Thus  Equation 
(6.8. 2.4)  can  be  written 
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v(t)  =  v(t)  -  v(o)  +  Iq  e^(t) 


(6. 8. 2. 5) 


It  should  be  emphasized  that  as  a  result  of  Option  6,  only  v(t)  is  actually 
calculated.  The  quantities  on  the  right  hand  side  of  Equation  (6.8. 2. 5)  are  not 
actually  known. 

When  Option  7  is  executed,  the  values  of  v(t)  on  the  DVSF  are  used  to 
calculate  Ig^t)  or  3(t),  the  derived  relative  displacements  found  from  v(t). 
Thus,  from  Equation  (6.8. 2. 5)  we  can  write 


d(t)  =  I^v(t)  =  I^v(t)  -  I%(o)  +  ihl  ^(t)  (6.8. 2.6) 

0  0  0  0  0 

=  d(t)  -  d(o)  and  v(o)  =  v(o)t.  Equation  6.8. 2.6  can  be  written 


Since  v(t) 


a(t)  =  d(t)  -  d(o)  -  v(o)t  +  ijlj  ^(t) 


(6. 8. 2. 7) 


d(t)  represents  true  displacement  at  time  t,  and  3(t)  approximate  relative 
displacement  at  time  t  found  from  v(t).  In  equations  (6.8. 2. 5)  and  (6. 8. 2. 7), 
the  actual  unknowns  are  v(t)  and  d(t),  respectively.  The  equations  can  thus  be 


re-wn  tten 


v(t)  =  v(  t)  -  [I|^  £^(t)  -  v(o) 


(6.8. 2.8) 


d(t)  =  d(t)  -  [IqIq  c^(t)  -  v(o)t  -  d(o) 


(6.8. 2. 9) 


The  quantities  in  square  brackets  are  then  correction  terms  that  must  be  applied 
to  extract  the  true  values  of  v(t)  and  d(t).  Let 


I'^(t)  =  7^(t)  -  v(o) 


:6.8.2.10) 
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and 


e^(t)  =  iJIq  e^(t)  -  v{o)t  -  d(o)  (6.8.2.11) 

Then  equations  (6.8. 2.3),  (6.8. 2.8)  and  (6.8. 2. 9)  can  be  written 

a(t)  =  a(t)  -  c^(t)  (6.8.2.12) 

v(t)  =  v(t)  -  e''(t)  (6.8.2.13) 

d(t)  =  d(t)  -  c^^lt)  (6.8.2.14) 

The  correction  procedure  thus  involves  determining  ^(t),  e'^(t),  and  e‘^(t)  so 
that  the  true  signals  can  be  extracted  fran  the  approximate  ones. 

Before  proceeding  further  a  simple,  but  most  common,  case  will  be 
considered.  Suppose  e^(t)  is  just  constant  in  time;  that  is,  e^(t)  =  where 
is  a  constant.  Then  Equation  (6.8.2.10)  shows  that 

e''(t)  =  ij,t  -  v{o)  (6.8.2.15) 

and  Equation  (6.8.2.11)  shows  that 

■“rl  12 

e  (t)=Yi|/t  -  v(o)t  -  d(o)  (6.8.2.16) 

~  V  ~  d 

Thus  the  e  (t)  is  just  a  polynomial  of  order  1,  a  linear  variation,  and  e  (t)  is 
a  polynomial  of  order  2,  a  quadratic  variation. 

For  generality  it  will  now  be  assumed  that  e'^(t)  can  be  represented 
adequately  as  a  polynomial  of  order  m: 

-H  -H  1  1 

e  (t)  =  )  (6.8.2.17) 

J=1  ^ 
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Once  the  coefficients  Cj  have  been  found.  Equation  (6.8.2.14)  can  be  used  to  find 
an  approximation  to  the  actual  displacements  d(t).  In  order  to  calculate  the 
coefficients  it  is  proposed  to  use  a  least  squares  procedure  as  outlined  in 
Appendix  C.  Now,  on  completion  of  Option  7,  a  DDSF  will  have  been  created 
containing  3^  values,  where  is  the  approximate  relative  displacement  at  the 
N'th  PN.  Let  a  least  squares  polynomial 


3(t) 


J-1 

J=1 


(6.8.2.18) 


be  fitted  to  N  data  points  starting  at  PN  P^  (for  t  =  0)  on  the  DDSF.  Thus  we 
want  to  fit  the  polynomial  in  the  least  squares  sense  to  the  data  points  (tp 

3n  „)  for  K  =  0toN-l.  It  should  be  recalled  that  the  t^  values  can  be 

calculated  from  Equation  (6.6. 2.1)  as 


'^P  +K  ^  (6.8.2.19) 

0 

where  h  is  the  reciprocal  of  the  digitizing  rate.  From  equations  (C3.5)  to 
(C3.7)  in  Appendix  C,  we  have  to  solve  the  system  of  equations 

m+1  _ 

,  A,,C',  =  Bj  (6.8.2.20) 

jLi  i 

for  I  =  1,  2,  ...  m+1 


where 


A 


IJ 


N-1 


K=0 


I+J-2 
^P  +K 


(6.8.2.21) 
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and 


=  =  I-l 

^I  ~  +K^P  +K 


(6.8.2.22) 


Since  the  (tp  dp  values  are  available  on  the  DDSF,  Ay  and  Bj  can  be 
calculated  and  Equation  (6.8.2.20)  solved  for  the  m+1  coefficients  Cj. 

The  next  step  is  to  derive  the  coefficients  £j  in  Equation  (6.8.2.17)  from 
the  coefficients  C,.  If  Equation  (6.8.2.14)  is  used  to  substitute  for  do  in 

Equation  (6.8.2.22)  and  then  Equation  (6.8.2.17)  is  used  to  substitute  for 

^0  ^ 

it  can  be  shown  after  a  straightforward  reduction  that 

m+l_  . 

I  Ay(Cj  -  Cj)  =  Bj  (6.8.2.23) 


for  I  =  1  to  m+1 


where 


N-1 

B,  =  H 

A  1/  _n 


l-l 


K=-0  Po-^^  Po"^ 


(6.8.2.24) 


Displacements  in  structural  response  are  often  oscillatory  about  a  zero 
displacement  baseline  as  shown  in  Figure  13.  Under  these  conditions,  the  moments 
Bj  will  oscillate  about  zero  when  considered  as  functions  of  N.  Although  all  the 
8j  will  not  in  general  be  zero  for  the  same  N,  it  will  usually  be  the  case  in 
practice  that  an  average  N  can  be  found  such  that  all  the  B.  are  sufficiently 

i  —  —  H 

close  to  zero.  It  will  then  be  true  from  Equation  (6.8.2.23)  that  C,  >=  C,  for 

I  csd  JO 

J  =  1  to  M  +  1.  e  (t)  can  then  be  evaluated  from  Equation  (6.8.2.17)  and  the 
correction  applied  to  all  points  d  on  the  DDSF  to  find  d(t). 

I 
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6.8.3  SPADE  Processing 

Figure  14(a)  shows  a  typical  acceleration  trace  from  an  accelerometer. 

Figure  14(b)  shows  a  plot  of  the  derived  relative  displacements  {5}.  It  appears 
that  the  displacements  in  this  case  are  superimposed  on  a  parabolic  type  of 
background  variation  as  anticipated  in  Equation  (6.8.2.16)  for  the  case  of 
e^(t)=  I})  (a  constant  error  offset  of  the  acceleration  trace). 

The  following  procedure  is  followed  by  SPADE  for  this  option  after  the  usual 
preliminary  questions  regarding  the  Series,  etc.,  have  been  answered. 

a)  SPADE  requests  a  time  (t^)  to  define  the  range  of  samples  on  the  DDSF 
to  be  used  in  the  least  squares  fit  for  the  determination  of  the 
coefficients  Cj  for  J  =  1  to  tiH-l  in  Equation  (6.8.2.18).  t^  should  be 
chosen  to  be  5  (or  greater  if  possible)  and  at  an  integral  number  of  cycles 
from  t  =  0. 

b)  SPADE  will  next  ask  whether  m  =  2  for  the  polynomial  to  be  fitted  is 
acceptable.  This  should  usually  be  answered  "yes"  unless  there  is  a  reason 
to  choose  a  higher  or  lower  order.  The  program  then  proceeds  to  do  the 
least  squares  fit  with  equations  (6.8.2.20)  to  (6.8.2.22).  On  completion  of 
the  fit,  the  coefficients  C,  are  used  as  the  Cj  in  Equation  (6.8.2.17). 
Equation  (6.8.2.14)  is  then  used  to  calculate  an  improved  approximation  to 
d(  t) . 

The  polynomial  (6.8.2.17)  is  next  differentiated  with  respect  to  time 

.  V 

to  give  t  (t).  If  the  usual  option  of  using  a  second  order  polynomial  in 
item  "b"  was  chosen,  (t)  will  then  be  a  polynomial  of  order  1  (a  linear 
function).  Equation  (6.8.2.13)  is  then  used  by  SPADE  to  correct  the 
velocities  on  DVSF.  The  polynomial  for  i:'^(t)  is  then  differentiated  again 
to  give  e^(t).  If  ^'^(t)  is  linear  then  c^(t)  will  just  be  a  constant. 
Equation  (6.8.2.12)  is  then  used  to  calculate  a(t). 
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Figure  15(a)  and  15(b)  show  the  results  for  the  case  in  Figure  14(b). 

In  this  case  ip  in  Equation  (6.8.2.16)  turned  out  to  be  only  about  2  G  [9]. 
This  amount  is  much  less  than  the  accuracy  to  be  expected  from  the 
accelerometers  but  leads  to  a  significant  improvement  in  the  results  as 
shown  in  Figure  15(b). 

If  a  graphics  terminal  is  being  used,  SPADE  will  optionally  plot  out 
the  graphs  shown  in  Figures  14(b)  and  15(b)  for  the  displacements.  The 
graphs  for  velocity  will  also  optionally  be  plotted  both  before  and  after 
correction. 

c)  On  examining  the  graphs  the  user  must  next  decide  whether  the  predicted 
initial  conditions  for  velocity  and  displacement  are  reasonable.  Referring 
to  equations  (6.8.2.10)  and  (6.8.2.11),  for  example,  it  is  seen  that  the 
least  squares  fit  to  c^(t)  will  have  produced  values  for  v(o)  and  d(o). 

SPADE  will  ask  whether  the  initial  conditions  are  satisfactory.  If 
not,  the  program  will  ask  whether  v(o)  is  to  be  specified.  If  "yes"  is 
given,  the  program  will  request  v(o),  apply  the  required  vertical  shift  of 
the  baseline  on  the  DVSF,  and  then  apply  a  correction  to  the  values  on  the 
DDSF  to  account  for  the  change  in  the  term  v(o)t  in  Equation  (6.8.2.11). 

The  program  will  next  ask  whether  d(o)  is  to  be  changed.  If  "yes"  is  given, 
the  program  will  request  d(o),  apply  the  required  vertical  shift  of  the 
baseline  on  the  DDSF,  and  then  apply  a  correction  to  the  values  on  the  DDSF 
to  account  for  the  change  in  d(o)  in  Equation  (6.8.2.11). 

In  conclusion,  the  correction  procedure  offered  in  Option  8  has  been  found 
to  be  considerably  effective  in  correcting  displacements  found  from  integrating 
acceleration  signals.  Although  the  discussion  above  assumed  elastic  response, 
the  present  version  of  SPADE  has  additional  procedures  that  allow  corrections  to 
be  applied  for  large  deformation  plastic  response. 
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6.9  Option  9:  FFT  and/or  Filtering  of  Transducer  Signals 

The  Fast  Fourier  Transform  Algorithm  (FFT)  originated  with  J.W.  Cooley  and 
J.W.  Tukey  [10].  SPADE  uses  a  Fortran  subroutine  presented  by  E.A.  Robinson  [11] 
to  calculate  the  Fourier  Transform  of  a  specified  transducer  signal.  A  sample 
run  of  SPADE  for  this  option  is  presented  in  Section  2.9  of  Reference  5. 
SPADE-PLOT  is  then  used  to  plot  the  transducer  time  signal  and  the  associated 
amplitude  ratio  against  frequency. 

6.10  Option  10:  Create  Averaged  Pressure  File  and  VAST  Load  File 

The  purpose  of  Option  10  is  to  create  a  load  file  for  the  program  VAST  [1] 
from  pressure  transducer  signals.  This  option  will  generate  a  pressure  loading 
for  each  element  of  a  finite  element  mesh.  A  sample  run  is  presented  in  Section 
2.10.2  of  Reference  5. 

6.11  Option  11:  Extract  Displacement,  Velocity,  and  Accewleration  from  VAST  file 

for  Plotting 

During  the  execution  of  the  VAST  finite  element  code,  quantities  such  as 
displacement,  velocity,  and  acceleration  will  be  copied  to  a  file  designated  as 
XXXX_T52,  where  XXXX  is  a  four  character  VAST  file  prefix.  This  file  will 
contain  kinematic  variables  for  all  degrees  of  freedom  at  every  node  and  thus 
could  be  quite  large. 

For  the  purposes  of  the  SPADE  program  it  is  convenient  to  introduce  the 
convention  of  a  Finite  Element  Transducer  Channel.  Tranducer  channels  were 
discussed  in  Section  3.4  and  depicted  in  Figure  2.  A  Finite  Element  Channel  is 
assigned  numbers  in  excess  of  500  and,  for  the  purposes  of  SPADE,  will  be  assumed 
connected  to  a  fictitous  transducer  which  produces  a  signal  that  corresponds  to 
the  finite  element  predictions.  Although  Finite  Element  Tranducers  can  be 
defined  independently  of  any  real  transducer,  for  correlation  studies  they  will 
usually  be  defined  to  be  spatially  coincident  with  some  real  transducer.  In  this 
case,  it  is  convenient,  although  not  necessary,  to  define  the  FE  Transducer  Code 
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to  be  the  same  as  the  real  one.  Thus  in  Figure  2,  for  example,  the  real 
transducer  A2  is  connected  to  channel  7.  The  corresponding  FE  Transducer  would 
have  the  same  Transducer  Code,  A2,  but  the  Channel  Number  would  be  set  to  507. 
Other  quantities,  such  as  the  Transducer  Coordinates,  would  then  be  the  same  for 
both  real  and  FE  transducers. 

The  purpose  of  Option  11  is  to  create  a  Finite  Element  Signal  Data  File 
(FESDF)  from  the  file  XXXX_T52  created  by  VAST.  The  FESDF's  have  the  same 
format  as  the  scaled  Signal  Data  Files  and  are  given  names  according  to  the  file 
naming  convention  in  Appendix  B2.  Once  the  FESDF  has  been  created,  it  can  be 
processed  and  plotted  using  SPADE  and  associated  plotting  programs.  The  FESDF  is 
essentially  considered  to  contain  a  signal  from  the  fictitous  Finite  Element 
Transducer. 


7.  DISCUSSION  OF  A  TYPICAL  RUN  OF  SPADE 


The  purpose  of  Section  7  is  to  summarize  briefly  the  steps  involved  in  using 
the  SPADE  program  for  processing  digitized  data  files  from  an  experiment. 
Reference  5  presents  considerably  more  detail  and  is  intended  to  guide  the  reader 
explicitly  through  a  series  of  sample  runs. 

The  first  step  in  processing  the  data  from  an  experiment  is  to  create  the 
Transducer  List  File  (TLF)  and  the  Transducer  Database  File  (TDF)  with  Option  1. 
Sections  B3  and  B4  in  Appendix  B  describe  the  format  of  these  files.  In 
particular.  Table  B4.1  gives  the  format  of  the  actual  data  stored  on  the  TDF. 
These  files  will  be  given  a  name  according  to  the  convention  described  in 
Section  B2 

Prior  to  running  Option  1,  the  COS  shown  in  Table  5.1  should  be  completed. 
Excuting  SPADE  will  then  allow  considerable  flexibility  in  the  way  data  can  be 
entered. 
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When  Option  1  has  been  completed,  the  next  step  is  to  determine  the  Signal 
Scale  Factors  from  the  Calibration  Pulses  with  Option  2.  To  do  this  for  a 
transducer,  SPADE  will  first  determine  the  CPDP  as  discussed  in  Section  6.2.3. 
This  will  then  be  followed  by  a  determination  of  the  calibration  pulse  height,  h, 
from  equation  6. 2.4. 2,  once  the  order  of  the  polynomial  for  the  background  has 
been  specified.  Equation  6. 2. 1.1  is  then  used  to  determine  E  from  the  CALIBRATE 
stored  on  the  TDF. 

When  all  the  Signal  Scale  Factors  have  been  determined.  Option  3  can  be  run 
for  each  transducer  to  obtain  the  Scaled  Signal  Data  Files.  This  will  involve 
correcting  the  signals  for  drift  and  then  multiplying  all  the  data  values  by  the 
appropriate  Signal  Scale  Factors. 

The  next  step  in  the  processing  should  be  to  synchronise  all  the  transducer 
signals  in  time  with  Option  4.  This  is  done  by  determining  a  point  on  each  SSDF 
that  corresponds  to  a  common  real  time  value.  This  common  time  value  is 
arbitrarily  associated  with  t  =  0.  In  order  to  make  the  choice  of  t  =  0 
physically  meaningful,  SPADE  offere  two  choices:  a)  the  average  time  of  arrival 
of  the  pressure  pulses  for  selected  pressure  transducers,  or  b)  the  first  time  of 
arri val . 

Option  5  should  next  be  executed  to  determine  the  pressure  pulse 
characteristics  such  as  the  positive  duration,  rise  time,  decay  constant,  and 
maximum  impulse.  Interactive  graphics  should  be  used  in  this  option  for  the  best 
resul ts. 

Option  6  should  next  be  executed  for  all  accelerometers  to  obtain  the 
relative  velocity  files.  Option  7  should  then  be  executed  similarly  to  obtain 
the  relative  displacement  files.  Finally,  Option  8  should  be  executed  to  obtain 
the  corrected  acceleration,  velocity,  and  displacement  files.  In  addition  to 
offering  several  types  of  correction  procedures  for  drift,  SPADE  also  allows,  if 
required,  nonzero  initial  conditions  to  be  applied  to  the  velocity  and 
displacement  files.  The  results  of  Options  6  and  7  are  thus  only  intermediate; 
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on  completion  of  Option  8  it  is  intended  that  the  absolute  velocities  and 
displacements  will  have  been  calculated  from  the  accelerations. 

Options  9  to  11  have  been  included  essentially  for  convenience  and  do  not 
necessarily  have  to  be  executed.  Detailed  examples  are  given  in  Reference  5. 


8.  CONCLUDING  REMARKS 


The  main  purpose  of  SPADE  is  to  give  a  facility  for  performing  basic  signal 
processing  on  the  large  volumes  of  data  created  during  blast  wave-structure 
interaction  experiments.  At  the  time  of  writing,  the  program  has  successfully 
been  used  by  five  term  students  to  process  a  considerable  amount  of  data.  The 
students  indicated  very  little  difficulty  in  using  the  code. 
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Figure  1 

BLOCK  DIAGRAM  OF  THE  DATA  AQUISITION,  PROCESSING, 
PLOTTING,  AND  ANALYSIS  PROCEDURES 
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Figure  2 

PICTORIAL  REPRESENTATION  OF  A  SIMPLE  SQUARE  PLATE 
INSTRUMENTED  WITH  A  PRESSURE  GAUGE  (P3),  A  STRAIN 
GAUGE  (S15,  AND  AN  ACCELEROMETER  (A2) 
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ILLUSTRATION  OF  A  CALIBRATION  PULSE  ON  A  CALIBRATION 
PULSE  DATA  FILE  WITHOUT  BACKGROUND  DRIFT 
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Figure  6 


ILLUSTRATION  OF  A  CALIBRATION  PULSE  WITH 
BACKGROUND  DRIFT  AND  NOISE 
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ILLUSTRATION  OF  A  TRANSDUCER  SIGNAL  STORED  ON  A 
SIGNAL  DATA  FILE  WITH  BOTH  DRIFT  AND  NOISE 
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ILLUSTRATION  OF  QUANTITIES  USED  FOR  THE  DETERMINATION 
OF  THE  START  TIME  (Tq)  FROM  PRESSURE  TRANSDUCER  SIGNALS 
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Figure  10 

ILLUSTRATIOIM  OF  PRESSURE  PULSE  CHARACTERISTICS 
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PLOT  OF  ACCELERATION  AGAINST  TIME  FOR  THE  TRANSDUCER  A  2 


Figure  14(b) 

PLOT  OF  DISPLACEMENT  AGAINST  TIME  FOR  THE  TRANSDUCER  A  2 


TYPICAL  PLOTS  OF  ACCELERATIOIM  AND  DERIVED  DISPLACEMENT 
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Figure  15(b) 

PLOT  OF  DISPLACEMENT  AGAINST  TIME  FOR  THE  TRANSDUCER  A 2 

TYPICAL  PLOTS  OF  CORRECTED  ACCELERATION  AND 
CORRECTED  DERIVED  DISPLACEMENT 
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APPENDIX  A 


The  EDIG  Data  Sheet  (EDS) 


A1  Introduction 


Two  sample  EDIG  data  sheets  are  shown  as  Tables  A1  and  A2.  The  header 
information  in  each  case  consists  of  the  following  information.  An  asterisk 
after  the  name  idicates  that  the  quantity  is  also  given  on  the  CDS  (Table  5.1). 

A2  Header  Information 


Charge  Size*.  See  Section  5. 

Height  of  Burst*.  This  is  referred  to  as  the  Standoff  on  the  CDS.  See 
Section  3.1  and  Fig.  2. 

Tape  Number.  This  is  an  identification  number  for  the  magnetic  tape  upon 
which  the  data  is  recorded. 

Date*.  This  is  the  date  on  which  the  experiment  was  conducted. 

Shot*.  See  Section  3.5, 

Ambient  Temp.  This  is  the  ambient  temperature  as  reported  by  the  DRES 
Meteorology  Section, 

Ambient  Press.  This  is  the  ambient  pressure  as  reported  by  the  DRES 
Meteorology  Section. 

Record  Speed.  This  indicates  the  speed  in  inches  per  second  (ips)  for  the 
magnetic  tape  when  recording  the  transducer  signals. 
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EXAMPLE  OF  AN  EDIG  DATA  SHEET  FOR  HONEYWELL  RECORDER 


EXAMPLE  OF  AN  EDIG  DATA  SHEET  FOR  AMPEX  RECORDER 


RESPONSE  OF  METAL  PLATES  TO  BLASE  WAVE  LOADING 
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EXAMPLE  OF  AN  EDIG  DATA  SHEET  FOR  AMPEX  RECORDER 
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Tape  Start.  This  is  the  number  of  feet  of  tape  that  has  been  allowed  to 
pass  the  Tape  Head  (Fig.  4)  prior  to  starting  the  tape  to  record  the  transducer 
signals  for  the  experiment.  The  number  of  feet  is  essentially  determined  by  the 
amount  of  information  already  on  the  tape  from  earlier  experiments.  The  IRIG  B 
time  may  also  be  given  for  the  point  at  which  the  tape  is  started  to  record  the 
experimental  data.  This  time  will  be  specified  in  the  same  format  as  for  the 
Event  (see  Section  5). 

Remarks.  This  is  for  additional  information. 

Record  Std.  This  is  an  IRIG  standard  of  recording.  WB  means  "Wide  Band". 
Several  different  recording  standards  exist. 

Event*.  The  number  of  feet  of  the  Tape  Head  (Fig.  4)  from  the  start  of  the 
tape  is  given  for  the  time  when  the  charge  is  exploded  (Det.  Zero).  The  IRIG  B 
time  at  Det.  Zero  is  also  given  (see  Section  5). 

Tape  Recorder.  This  will  indicate  either  an  Ampex  or  Honeywell  recorder  as 
discussed  in  Section  3.4. 

Tape  Stop.  This  is  the  number  of  feet  of  tape  that  has  passed  the  Tape  Head 
(Fig.  4)  when  the  tape  is  stopped  after  the  experiment.  This  amount  can  then  be 
used  to  position  the  Tape  Head  for  a  future  experiment.  The  IRIG  B  time  may  also 
be  given  corresponding  to  the  time  at  which  the  tape  is  stopped. 

A3  Columns  on  the  EDS 

A3.1  Introduction 


There  are  a  total  of  15  columns  on  the  EDS.  Columns  1  to  3  contain  gauge 
data.  Columns  4  to  6  contain  transducer  location  information,  Columns  7,  8,  14 
and  15  contain  additional  information,  anc  Columns  9  to  13  contain  signal 
conditioning  information. 
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A3. 2  Gauge  Data  (Columns  1  to  3) 

Column  1  contains  the  Transducer  Codes  as  discussed  in  Section  3.2. 

Column  2  contains  the  Transducer  Model  and  Serial  Numbers.  Thus  PI  has 
Model  Number  113A21,  Serial  Number  2686.  Also  P2  to  P4  have  the  same  model 
number  (113A21)  and  the  indicated  Serial  Numbers. 

Column  3  contains  certain  gauge  constants.  These  are  discussed  fully  in 
Reference  5  along  with  signal  conditioning. 

A3. 3  Transducer  Location  Information  (Columns  4  to  6) 


Columns  4  to  6  contain  transducer  location  data.  Column  4  gives  the  surface 
on  which  the  transducer  is  mounted  and  is  the  value  listed  in  Column  6  of  the  CDS 
(Table  5.1).  Column  5  contains  the  Transducer  Coordinates  (XC,  YC)  discussed  in 
Section  3.3  and  listed  in  Columns  3  to  4  of  the  CDS.  For  strain  gauges  only, 
Column  6  contains  the  coordinate  ANGLE  given  in  Column  5  of  the  CDS. 

A3. 4  Additional  Information  (Columns  7,  8,  14  and  15) 

Column  7  contains  just  the  estimated  maximum  levels  of  stimulus  to  which 
each  gauge  is  expected  to  be  subjected.  Thus  PI  to  P3  are  expected  to  be 
subjected  to  pressures  of  10  PSI.  A2  is  expected  to  be  subjected  to  an 
acceleration  of  970G.  S7  is  expected  to  read  a  maximum  compressive  strain  of 
80ue,  where  ye  denotes  microstain.  The  values  in  Column  7  are  used  by  EDIG  to 
set  the  gain  settings  for  signal  conditioning  as  will  be  explained  in  Section 
A3. 5. 


Column  8  contains  cable  identification  numbers  used  to  physically  tag  each 
of  the  cables  attached  to  transducers.  The  tags  correspond  to  terminal  numbers 
which  allow  EDIG  to  identify  the  transducers  attached  to  each  cable.  It  should 
be  mentioned  that  these  numbers  are  not  the  Channel  Numbers  referred  to  in 
Section  3.4. 
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Column  14  contains  the  track  numbers  on  the  magnetic  tape  used  for  the 
transducer  signals  (see  Section  4.2  and  Fig.  4).  Table  A1  for  the  H-101  tape 
recorder  indicates  that  the  signals  from  transducers  Pi  to  P4  were  stored  on 
tracks  1  to  4,  the  signals  from  transducers  A1  to  A6  were  stored  on  tracks  5  to 
11,  and  the  signals  from  transducers  S7  and  S8  were  stored  on  tracks  12  and  13, 
respectively.  Track  14  was  used  to  store  the  IRI6  signal.  Similarly,  Table  A2 
shows  the  track  numbers  assigned  to  transducers  S9  to  S22  on  the  Ampex  tape. 

Note  that  on  both  the  Honeywell  and  Ampex  tape,  track  13  was  used  to  record  an 
actual  transducer  signal  rather  than  the  Det.  Zero  voltage  pulse  referred  to  in 
Section  4.2.  This  was  done  in  this  case  to  allow  26  rather  than  24  transducer 
signals  to  be  recorded.  All  the  transducer  signals  were  frequency  modulated 
(F.M.)  to  reduce  the  background  noise  effects. 

Column  15  is  used  to  record  additional  information. 

Additional  information  on  the  EDIG  data  sheet  can  be  found  in  Appendix  A  of 
Reference  5. 
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APPENDIX  B 

SPADE  Data  and  Database  Files 


B1  Introduction 


Several  types  of  files  have  been  referred  to  throughout  this  report.  These 
are  the  Transducer  Database  File  (TDF),  Transducer  List  File  (TLF),  Calibration 
Pulse  Data  File  (CPDF),  Signal  Data  File  (SDF),  Scaled  Signal  Data  File  (SSDF), 
Derived  Velocity  Signal  File  (DVSF),  and  Derived  Displacement  Signal  File  (DDSF). 
In  the  next  section,  the  file  naming  convention  for  each  will  be  discussed.  The 
following  sections  will  then  describe  the  contents  and  data  format. 

B2  File  Naming  Convention 

For  the  purpose  of  describing  the  file  naming  convention  the  following 
notation  will  be  used. 


Table  B2.1 

Symbols  Used  to  Describe  the  File  Naming  Convention 


Symbol 


Meani ng 


A-B- 

www 

XX 

YY 

11 


Cnn 


User's  prefix. 

A  Channel  Number  of  up  to  three  digits. 

A  Series  Number  of  up  to  two  digits. 

A  Shot  Number  of  up  to  two  digits. 

A  Panel  Number  of  up  to  two  digits. 

A  code  which  is  A  for  acceleration,  V  for  velocity,  U 
for  displacement,  S  for  strain,  P  for  pressure,  and  I 
for  impulse. 

A  Transducer  Code  (see  Section  3.2)  for  a  transducer 
which  measures  a  stimulus  of  type  C,  where  C  is  a 
single  character  assigned  according  to  the  same 
convention  as  T.  nn  is  an  integer  of  up  to  two  digits. 
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Table  B2.2  below  gives  the  file  naming  convention  for  each  type  of  file. 
The  overscored  parts  refer  to  the  notation  described  in  Table  B2.1. 


Table  B2.2 


File  Type 


File  Name  Convention 


A-B-TRANLIST-SEXXSHYY 


A-B-TDATA-SEXXSHYY 


CPDF  recorded  on 
H-101  (Section  3.4 
and  Fig.  2) 


H  CCHWWWSHYY  XX 


CPDF  recorded  on 
A-PR2200  (Section 
3.4  and  Fig.  2) 


CCHWWWSHYY  XX 


SDF  recorded  on 
H-101  (Section  3.4 
and  Fig.  2) 


H  CHWWWSHYY  XX 


SDF  recorded  on 
A-PR2200  (Section 
3.4  and  Fig.  2) 


CHWWWSHYY  XX 


A-B-SEXXPZZSYYCnn^T 
(C  and  T  must  desTgnate  the 
same  type  in  this  case) 

A-B-SE^PZZSYYCnn_V 
(C  will  denote  the  type  of 
transducer  from  which  the 
velocity  was  derived  and 
will  usually  be  A) 


A-B-SEXXPZZSYYCnn_U 
(C  will  denote  the  type  of 
transducer  from  which  the 
displacement  was  derived 
and  will  usually  be  A). 
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B3  Transducer  List  Files 


The  TLF  are  created  during  the  execution  of  Option  1  in  SPADE.  The  purpose 
of  the  TLF  is  to  summarize  the  transducers  that  are  on  the  TDF  and  is  used  for 
display  purposes  only.  The  first  few  lines  of  the  file  contain  header 
information.  After  this  a  listing  of  the  Channel  Numbers,  Transducer  Codes, 
Surface  Codes,  and  Transducer  Coordinates  is  given. 

B4  Transducer  Database  Files 


The  TDF  will  contain  several  lines  of  data.  The  file  is  formatted  and  can 
thus  be  altered  with  the  CP6  editor  if  required.  One  of  the  lines  in  the  file 
will  contain  three  asterisks  "***''  in  Columns  1  to  3.  All  lines  prior  to  this  in 
the  file  will  be  treated  as  comment  lines  by  SPADE.  Thus  additional  comment 
lines  can  be  placed  in  this  region  with  the  CP6  editor  or  otherwise  if  required. 
All  lines  after  the  line  containing  the  asterisks  are  actual  data  lines. 

Each  line  of  data  in  the  TDF  contains  the  quantities  listed  in  the  following 
table.  The  "Terminology"  and  "Nomenclature"  given  at  the  start  of  this  manual 
can  be  referred  to  for  clarification  if  required.  The  "Terminology"  also  gives 
section  references. 
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Table  B4.1. 


Quantities  Stored  on  the  TDF 


Quantity 

Columns 

Quanti ty 

Columns 

CHANNEL 

1 

- 

5 

CPDP  P^ 

b 

92 

-  95 

TRANSDUCER  CODE 

7 

- 

9 

SDP  Pj 

97 

-  100 

CALIBRATE 

11 

- 

22 

SDP  P^ 

102 

-  105 

SIGNAL  SCALE  FACTOR 

24 

- 

35 

SDP  P3 

107 

-  no 

DRATE 

37 

- 

48 

SDP  P^ 

112 

-  115 

XC 

50 

- 

53 

SDP  P^ 

117 

-  120 

YC 

55 

- 

58 

SDP  P. 

b 

122 

-  125 

ANGLE 

60 

63 

P  (Point  Number 

127 

-  130 

for  t  =  0) 

SURFACE 

65 

1 

CPDP  P^ 

67 

- 

70 

Decay  Constant 

132 

-  143 

CPDP 

72 

- 

75 

Maximum  Impulse 

145 

-  156 

CPDP  P3 

77 

- 

Digitizing  Start  Time 

CPDP  P, 

82 

85 

Y 

157 

-  160 

D 

161 

-  163 

CPDP  P5 

87 

- 

90 

H 

164 

-  165 

M 

166 

-  167 

S 

168 

-  169 

MS 

170 

-  172 

MICRSEC 

173 

-  175 

i 

Pf,  (for 

177 

-  188 

curve  fit) 

In  Table  B4.1,  the  Channel,  Transducer  Code,  Calibrate,  Drate,  XC,  YC, 
Angle,  Surface,  and  Digitizing  Start  Time  are  placed  in  the  TDF  during  the 
execution  of  Option  1  as  described  in  Section  6.1.  An  exception  occurs  should  a 
CPDF  not  be  available.  In  this  case  the  Signal  Scale  Factor  is  determined  in 
Option  1  on  the  basis  of  the  Calibration  Constant  and  Gain  supplied  on  the  CDS. 
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The  Signal  Scale  Factor  and  CROP  (the  to  P^)  are  determined  in  Option  2 
except  as  noted  above  when  a  CPDF  is  not  available.  The  SDP  P^,  P^,  P^  and  P^ 
are  determined  in  Option  3.  P^  is  determined  in  Option  4.  The  Decay  Constant 
and  Maximum  Impulse  are  determined  in  Option  5,  as  well  as  SDP  Pj  and  P^. 

B5  CPDF  and  SDF 


The  CPDF  and  SDF  contain  data  samples  from  the  analog  tape  (see  Section  4 

and  Fig.  4).  The  format  for  each  line  is  1018.  Typically  there  will  be  say  5000 

samples  on  the  file  (and  thus  500  lines).  As  explained  in  Section  4.3,  with  a 

GAIN  of  1,  a  1  Volt  signal  on  the  analog  tape  will  produce  the  number  204.8 

during  the  digitizing  process.  This  would  be  stored  on  the  CPDF  or  SDF  as  205  in 
the  format  18. 

86  SSDF,  DVSF  and  DDSF 


The  SSDF.  DVSF  and  DDSF  are  created  by  SPADE  during  processing.  The  format 
for  each  line  is  10E12.5.  In  effect,  the  SDF  data  stored  as  1018  is  replaced  by 
the  processed  data  stored  as  10E12.5 
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APPENDIX  C 


Theory  of  Least  Squares  Polynomial  Curve  Fits 


The  requirement  is  to  fit  a  polynomial  of  the  form 


y  =  I  C.x 


(C3.1) 


to  a  set  of  N  data  points  (x^,  yj^)  where  R  goes  from  1  to  N.  m  is  the  order  of 
the  polynomial  and  the  C^  are  the  m+1  coefficients.  The  following  condition 
hold; 


Condi tion 


Consequence 


1.  m  >  N-1 


A  unique  polynomial  cannot  be  fitted  to 
the  data  points. 


2.  m  =  N-1 

and  no  two  data  points 
have  the  same  x  value. 


A  unique  polynomial  can  be  fitted  that 
will  pass  through  all  the  data  points. 


3.  m  N-1 


There  is  no  polynomial  of  order  m  that 
will  pass  through  all  the  data  points. 


For  condition  3  it  is  possible  to  determine  a  polynomial  of  order  m  that 
will  approximate  the  values  of  the  data  points.  Using  the  theory  of  least 
squares,  this  is  done  by  minimizing  the  sum,  S,  of  the  square  of  the  deviations 
of  the  polynomial  from  the  data  points.  S  is  given  by 


N  m+1  ,  ,  2 

S  =  nyp  -  )  C, 

R  =  1  J  =  1  ^ 

It  is  then  necessary  to  minimize  S  by  varying  the  coefficients  Cj  of  the 
polynomial.  The  condition  for  a  minimum  is 


(C3.2 
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represents  a  column  matrix  containing  the  m+1  unknown  coefficients  Cj;  namely, 


(C3.9) 


[B}  also  represents  a  column  matrix  containing  the  m+1  values  of 


(C3.10) 


Provided  m+1  >  N  and  not  less  than  m+1  of  the  data  points  have  distinct  x 
coordinates,  the  system  of  equations  C3.8  can  be  solved  for  the  unknown 
coefficients  |C}.  In  the  special  case  m+l-N  (with  all  the  data  points  having 
unique  x  coordinates)  the  solution  of  the  system  of  equations  will  provide  the 
coefficients  |C|  that  define  the  unique  polynomial  that  passes  through  all  the 
data  points  according  to  condition  2  stated  above. 

The  system  C3.8  is,  in  general,  badly  conditioned  but  provided  m  <  5  and 
double  precision  is  used  in  the  calculations  the  accuracy  will  normally  be 
adequate.  SPADE  automatically  limits  m  to  the  range  stated.  Values  of  zero  for 
m  are  permisible  to  SPADE  also.  In  this  case  the  single  coeffcient  Cj  will  be 
the  average  of  the  yj  values. 
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